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ABSTRACT 


The  Sensor  Hosting  Autonomous  Remote  Craft  (SHARC),  also  known  as  Wave  Glider,  is 
an  autonomous  oeean  vehiele  powered  by  wave  motion.  This  slow-moving  platform 
makes  long-term  deployments  and  environmental  data  collection  feasible,  especially  in 
data  sparse  regions  or  hazardous  environments.  The  standard  SHARC  hosts  a 
meteorological  station  (Airmar  PB200)  that  samples  air  pressure,  temperature,  wind 
speed  and  wind  direction  at  1.12  m.  The  SHARC  automatically  transmits  a  10-minute 
averaged  data  suite  through  an  Iridium  satellite  link.  In  an  effort  to  evaluate  the  SHARC 
default  Airmar  sensors  and  seek  optimal  sensors  for  air-sea  interaction  studies,  NFS  has 
developed  an  independent  package  of  meteorological  sensors,  the  “NFS  Met,”  for  use  on 
the  SHARC.  NPS  Met  measures  pressure,  air  temperature,  wind,  SST,  and  relative 
humidity.  This  SHARC  payload  package  was  deployed  three  times  in  the  Monterey  Bay, 
along  with  a  collocated  drifting  buoy  (Marine  Air-Sea  Flux  buoy,  or  MASFlux)  with 
proven  flux,  mean,  wave,  and  SST  measurement  for  comparison  and  validation.  This 
thesis  will  present  analyses  of  data  from  the  new  mast  and  Airmar  as  compared  to  known, 
quality  measurements  from  NPS  MASFlux  and  NDBC  buoy.  Surface  fluxes,  evaporation 
duct  heights  and  strength  are  derived  from  the  SHARC  measurements  using  the  COARE 
algorithm. 
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I.  INTRODUCTION 


A.  IMPORTANCE  OF  NEAR  SURFACE 

Accurately  quantifying  the  eharaeteristies  of  the  near  surfaee  environment  is 
eritieal  to  understanding  and  modeling  Earth’s  atmosphere.  In  partieular,  studying  the 
near  surfaee  marine  environment  is  espeeially  important  beeause  oeeans  eover  over  70 
pereent  of  the  surface  of  the  Earth.  The  exehanges  of  momentum,  heat,  moisture,  and 
turbulenee  among  the  upper  oeean,  oeean  surfaee,  and  lower  atmosphere  play  an  integral 
role  in  determining  the  dynamics  of  both  the  oeean  and  the  atmosphere.  Better  theoretieal 
representation  and  experimental  examination  of  the  interplay  between  the  oeean  and 
atmosphere  are  needed,  whieh  ealls  for  a  major  systematie  and  well-organized 
observation  program  of  the  oeean  and  the  strueture  of  the  atmospherie  boundary  layer 
(Kraus  and  Businger  1994).  Obtaining  researeh  grade,  high  quality  near  oeean  surfaee 
measurements,  however,  poses  many  ehallenges.  A  long  enduranee,  mobile  platform 
eapable  of  observing  the  near  surfaee  environment  would  benefit  many  scientifie  fields. 

Collecting  in  situ  measurements  of  key  variables  that  eharaeterize  the  nonlinear 
air-sea  interaetion  proeess,  or  an  aeeurate  representation  of  sueh  proeesses,  remains  one 
of  the  greatest  ehallenges  to  environmental  foreeast  model  improvement.  Eor  the 
exehanges  of  mass  and  energy  aeross  the  air-sea  interfaee,  bulk  aerodynamie  surfaee  flux 
parameterization  sehemes  are  needed,  which  are  largely  built  upon  the  Monin-Obukhov 
similarity  theory  (MOST)  developed  in  1954  (Eiu  et  al.  1979).  Eairall  et  al.  (1996b,  2003) 
further  refined  the  bulk  surface  flux  seheme  using  data  eolleeted  during  the  Tropical 
Ocean  Global  Atmosphere  (TOGA)  Coupled  Oeean-Atmosphere  Response  Experiment 
(COARE)  that  resulted  in  what  is  known  as  the  COARE  bulk  flux  algorithm.  COARE  is 
now  the  most  widely  used  flux  algorithm  in  most  of  the  mesoseale  and  global  seale 
foreeast  models.  Eimitations  to  the  COARE  algorithm  do,  however,  exist,  partly  due  to 
the  inherent  assumptions  of  MOST  (e.g.,  Andreas  et  al.  2014).  Thus,  the  need  to  eolleet 
high-quality  data  at  the  air-sea  interfaee  persists  in  order  to  further  improve  surfaee  flux 
parameterization  algorithms  and  henee  foreeast  models. 
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Near-surface  measurements  are  also  important  for  developing  and  verifying 
coupled  air-sea  forecast  models.  Coupled  air-sea  models  have  produced  improved  results 
in  certain  conditions  and  demonstrate  the  potential  for  greater  use  (Powers  and  Stoelinga 
1999).  Coupled  models  are  also  critical  for  studying  climate  variability  because  the 
accurate  parameterization  of  air-sea  momentum  and  heat  flux  is  critical  for  realistic 
simulation  of  the  atmosphere-ocean  system  (Brown  et  al.  2013).  Near  surface  data  sets 
and  observations  are  thus  indispensable  to  fill  in  the  need  of  developing  and  improving 
coupled  forecast  models. 

Near  surface  measurements  also  provide  critical  validation  of  satellite  data  that 
provide  broad  global  coverage  of  measurements,  particularly  in  otherwise  data  sparse 
regions.  Developing  and  refining  remote  sensing  techniques,  however,  require  in  situ 
observations.  Due  to  a  lack  of  open  ocean  near  surface  measurements,  many  satellite 
validations  over  the  ocean  rely  on  climatology  or  bulk  algorithms  to  convert  ship  or  land- 
based  measurements  to  near  surface  quantities.  These  data  sets  are  not  adequate  to  fully 
validate  satellite  retrieval  methods  over  the  open  ocean,  thus,  more  in  situ  observations 
are  necessary  (Waliser  et  al.  1999). 

Similarly,  electromagnetic  (EM)  propagation  models  also  depend  on  robust  in  situ 
measurements  and/or  highly  accurate  forecasts.  Surface  moisture  and  temperature 
gradients  affect  the  height  of  the  evaporative  duct  that  can  act  as  a  wave  guide  for  high 
frequencies  (Edson  et  al.  1999).  EM  propagation  is  very  sensitive  to  the  vertical  variation 
in  the  atmospheric  boundary  layer  (ABE)  and  would  also  benefit  from  improvements  of 
model  parameterizations  with  additional  data  for  both  research  and  operational  purposes. 

B,  RESEARCH  OBJECTIVES 

A  network  of  near  surface,  ocean  observing  platforms  that  are  capable  of 
providing  high  quality  measurements  would  aid  in  furthering  all  of  the  studies  described 
above.  In  2012,  NFS  received  two  wave  gliders,  referred  to  by  the  Navy  as  “Sensor 
Housing  Autonomous  Remote  Crafts”  (SHARC).  The  goal  of  this  thesis  is  to  improve  the 
quality  and  coverage  of  near  surface  observations  available  for  Naval  applications  and 
surface  flux  research.  Specifically,  this  thesis  will  present  an  evaluation  of  the  SHARC  as 
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a  platform  for  gathering  meteorological  and  oceanographic  (METOC)  data  in  general, 
and  air-sea  interface  flux  data  in  particular.  Compared  to  other  near  shore  platforms,  the 
SHARC  adds  longevity,  mobility,  and  the  ability  to  “drive”  to  areas  of  interest. 
Measurements  from  the  standard  SHARC  weather  station  sensors,  as  well  as  a  newly 
developed  suite  of  sensors,  will  be  compared  to  those  from  research  grade  sensors  on  the 
NFS  Marine  Air-Sea  Flux  (MASFlux)  buoy.  Three  co-deployments  of  the  SHARC  and 
the  MASFlux  buoy  in  Monterey  Bay,  and  static  tests  of  the  SHARC  standard  instruments 
and  NFS  sensors,  were  conducted  to  provide  the  datasets  for  comparison.  These  datasets 
will  also  be  used  for  surface  flux  parameterization  analysis.  Additionally,  other  SHARC 
users  have  provided  measurements  collected  by  the  standard  SHARC  weather  sensors. 
These  groups  have  presented  initial  findings  in  Wang  et  al.  (2012)  and  Griffith  et  al. 
(2012).  This  thesis  will  provide  a  more  detailed  examination  of  their  datasets. 

This  thesis  focuses  on  the  testing  and  improvement  of  the  SHARC  and  its 
environmental  sensors.  The  ultimate  goal  of  this  thesis,  however,  is  the  utilization  of 
SHARCs  to  provide  high  quality  near  surface  data  in  sparse,  hostile,  or  severe  weather 
environments.  The  resulting  data  can  then  be  used  in  applications  as  varied  and  important 
as  protecting  Naval  assets  and  personnel  from  adverse  weather,  to  achieving  decision 
superiority  in  the  electromagnetic  spectrum. 

C.  NAVAL  APPLICATIONS 

The  near  surface  environment  plays  a  critical  role  in  determining  EM  propagation, 
improving  coupled  model  and  severe  weather  forecasting,  and  in  predicting  wave 
dynamics  and  forecasts.  As  the  Navy’s  focus  shifts  from  mainly  blue  water  operations  to 
also  include  more  littoral  regions,  the  complexity  of  the  environment,  forecasting 
requirements,  and  decision  aids  also  increase.  Rear  Admirals  White  and  Filipowski  said 
in  a  recent  Proceedings  article,  “[cjollecting  the  raw  data,  processing  and  analyzing  the 
data,  and  predicting  how  that  data  will  change...  assist  operational  commanders  in 
making  faster  and  better  decisions  than  the  adversary”  (2014).  According  to  a  1999 
Office  of  Naval  Research  (ONR)  report,  the  Navy  requires  increasingly  accurate  and 
detailed  weather,  wave,  and  ocean  forecasts  (Edson  et  al.).  There  is  also  mounting 
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interest  in  eoupled  models  and  foreeasts.  Near  surfaee  data  are  key  to  improving  surfaee 
flux  parameterization  and  understanding  the  meehanisms  that  eouple  the  oeean  and 
atmosphere.  Additionally,  real-time  near  surfaee  observations  are  essential  for  model 
development  and  validation. 

1.  Electromagnetic  Propagation 

The  most  important  and  immediate  applieations  of  near  surfaee  measurements  is 
for  EM  propagation  studies  and  predietions.  The  Navy  uses  the  EM  speetrum  for 
eommunieation,  radar  deteetion  and  avoidanee,  and  eleetronie  warfare.  As  weapon 
systems  and  platforms  become  more  complex,  so  does  their  dependence  on  EM 
propagation  and  control  of  the  EM  domain.  EM  performance  predictions  are  therefore  a 
Navy  requirement.  As  Admiral  Jonathan  Greenert,  Chief  of  Naval  Operations,  said. 

The  electromagnetic  (EM)  spectrum  is  an  unseen  but  integral  part  of  our 
daily  personal  lives.  And,  it  is  essential  to  our  military  operations,  bailing 
to  use  it  effectively  WILE  be  the  difference  between  victory  and  defeat. 

The  EM  environment  presents  us  with  challenges  and  opportunities  in  the 
21st  century.  It  is  an  area  we  can  use  to  gain  an  advantage  over  our 
adversaries.  To  command  it,  we  need  the  ability  to  monitor  and  be  aware 
of  the  environment,  manage  our  emissions,  discretely  communicate,  find, 
track  and  defeat  threats,  and  conduct  attacks  as  needed.  (2012) 

As  described  above,  EM  propagation  is  very  sensitive  to  even  small  changes  in 
the  ABE  temperature  and  moisture  gradients.  Near  surface  measurements  are  critical  to 
determining  and  accurately  forecasting  these  gradients  and  the  structure  of  the  ABE. 
Advanced  Refractive  Effects  Prediction  System  (AREPS)  and  Target  Acquisition 
Weapons  Software  (TAWS),  in  particular,  require  accurate  ABE  profile  inputs  in  order  to 
effectively  predict  EM  propagation,  radar  ranges,  and  weapon  sensor  effectiveness.  In 
fact,  Konstanzer  (1994)  found  that  a  two  meter  change  in  the  evaporation  height  duct 
(EDH)  could  cause  a  50  percent  change  in  the  SPY-1  radar  detection  range.  Both  in  situ 
near  surface  measurements,  and  improved  parameterization  equations  will  improve  the 
Navy’s  EM  performance  prediction  capability. 
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II.  BACKGROUND 


The  interaction  between  the  atmosphere  and  the  ocean  at  the  air-sea  interface  is 
dynamic,  complex,  and  non-linear.  This  interaction,  and  the  associated  exchange  of 
momentum,  energy,  and  mass  across  the  air-sea  boundary,  determines  the  structure  of  the 
lower  atmosphere  and  upper  ocean.  These  processes,  while  critical  to  many  scientific 
fields  as  described  in  Chapter  I,  remain  difficult  to  measure  or  study  directly. 

Surface  flux  parameterizations  seek  to  model  these  air-ocean  exchanges.  A  further 
extension  of  the  surface  flux  parameterization  is  to  calculate  near  surface  profiles  when 
in-situ  measurements  of  fluxes  and  profiles  with  sufficient  vertical  resolution  are  not 
available.  The  state-of-art  COARE  algorithm,  validated  by  a  large  amount  of  over-ocean 
datasets,  has  been  broadly  used  by  nearly  all  mesoscale  forecast  models  as  well  as  models 
at  smaller  or  larger  scales.  Combining  with  MOST,  it  is  also  used  to  derive  surface  layer 
profiles  for  applications  such  as  EM  propagation  and  ducting  over  the  ocean.  In  this 
study,  deduction  of  surface  flux  and  surface  layer  vertical  profiles  from  the  SHARC 
measurements  will  also  continue  to  utilize  MOST  and  COARE. 

A.  PHYSICAL  PROCESSES  AT  THE  AIR-SEA  INTERFACE 

The  ocean-atmosphere  exchange  consists  of  many  dynamic  and  interwoven 
processes,  as  shown  in  Figure  1.  At  the  heart  of  this  interaction  are  the  turbulent 
exchanges  of  momentum,  sensible  and  latent  heat,  which  drive  turbulent  mixing  in  the 
atmospheric  boundary  layer.  Such  turbulent  mixing,  combined  with  entrainment  at  the 
top  of  the  boundary  layer,  modifies  the  near-surface  wind,  temperature,  and  specific 
humidity,  and  hence  feedback  to  air-sea  exchange  processes.  The  air-sea  momentum  and 
energy  transfer  is  thus  a  coupled  process.  On  the  upper  ocean  side,  wind-driven  surface 
stress,  sensible  and  latent  heat  fluxes,  and  atmospheric  radiation  result  in  upper  ocean 
turbulent  mixing  that  modifies  the  sea  surface  temperature.  Wind  generated  waves  and 
swells  are  also  part  of  the  equation  in  these  coupled  processes  to  modify  the  surface 
roughness.  The  wave  breaking  creates  aerosols  and  sea  spray  that  affect  the  turbulent  heat 
transfer  and  the  absorption  and  reflection  of  radiation.  Changes  in  radiation  and  heating 
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also  affect  precipitation,  which,  in  turn,  may  modify  the  wind  and  thermodynamics  of  the 
subcloud  layer  through  evaporative  cooling.  On  the  larger  scale,  these  air-sea  processes 
drive  open  ocean  circulation  and  determine  the  mixed  layer  depth.  Within  the 
atmosphere,  these  processes  determine  the  structure  of  the  lower  atmosphere. 
Understanding  near  surface  physical  processes  is  therefore  critical  to  modeling  and 
predicting  the  coupled  air-ocean  environment. 


Figure  1.  A  schematic  of  the  air-sea  interface  processes  (from  Woods  Flole 

Oceanographic  Institution  [WHOI]  2011) 

B,  SURFACE  FLUX 

The  energy,  mass,  and  momentum  exchange  between  the  atmosphere  and  the 
ocean  physically  couple  the  two  media.  The  amount  of  exchange  of  a  quantity  across  a 
unit  area  is  called,  “flux.”  Flux  occurs  as  a  result  of  flow  perpendicular  to  and  through  the 
interface.  For  example,  the  amount  of  flux  (F^)  for  a  variable  x  through  a  horizontal 
plane  by  the  vertical  wind,  w,  could  be  shown  by 


F^^wx 


(1) 
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The  flow  and  the  entity  being  examined  then  can  be  broken  into  mean  (denoted  by  the 
overbar)  and  perturbation  (denoted  by  the  prime)  portions  using  Reynolds  decomposition. 


=  wx  =  (w  +  w')»(  X  +  x')  =  wx  +  w' x' 

The  first  term  on  the  right-hand-side  of  Eqn.  (2)  is  referred  to  as  the  advective  flux 
caused  by  the  mean  flow  of  a  fluid,  while  the  second  term  is  a  result  of  turbulent  or 
perturbation  flow  and  is  referred  to  as  the  eddy  flux  or  turbulent  flux.  For  the  horizontal 
plane,  the  mean  vertical  velocity  is  generally  small,  thus  reducing  the  total  flux  to 
turbulence  flux  only. 


One  can  then  use  the  eddy  correlation  method  to  calculate  vertical  turbulent  flux 
of  a  quantity  as  the  covariance  of  the  perturbations  in  vertical  wind  velocity  and  the 
perturbations  in  the  desired  entity  at  a  point  (Readings  et  al.  1969).  With  sufficiently  fast 
sampling  instruments,  this  correlation  method  allows  for  direct  turbulence  flux 
measurements. 


The  most  important  fluxes  between  air  and  sea  at  the  interface  are  the  vertical 
exchange  of  momentum  (r),  sensible  heat  and  latent  heat  {Hi).  These  fluxes  can  now 
be  represented  by  scaling  parameters  defined  by  the  following  equations: 

T  =  p{w'u'  +  w'v'  )'^^  = 

//,  =  pCp  w' O'  =  -pCpUM  (4) 

=  pL^w'q'  =  -pL^u.q, 


where  p  is  the  dry  air  density,  u’,  v’,  w’  are  the  horizontal  and  vertical  wind  components; 
Cp  is  the  specific  heat  capacity  of  air  at  constant  pressure;  is  the  latent  heat  of 

vaporization;  6  is  air  potential  temperature;  the  prime  indicates  perturbations;  the  overbar 
indicates  the  time  average;  and  the  ‘  *  ’  subscript  indicates  the  respective  quantity  surface 
layer  turbulence  perturbation  scale. 

Direct  measurements  of  fluxes  and  their  governing  turbulent  processes  are 
difficult  to  obtain  operationally.  These  exchanges  happen  on  small  scales  and  require 
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high-rate  sampling  so  that  the  small  seale  turbulent  perturbations  ean  be  resolved. 
Therefore,  ealeulating  flux  by  using  more  easily  obtained  mean  profile  measurements  is 
eritieal  to  modeling  the  environment  and  the  numerous  applications  described  in  Chapter 
I.  This  concept  is  the  basis  for  surface  flux  parameterization. 


C.  MONIN-OBUKHOV  SIMILARITY  THEORY  (MOST) 

In  the  lowest  10  percent  of  the  atmospheric  boundary  layer,  Coriolis  is  negligible 
and  flux  remains  nearly  constant  with  height  (Businger  et  al.  1971).  This  layer  is  called 
the  surface  layer  and  is  approximately  50-1 00m  in  height.  In  this  layer,  the  primary 
sources  for  generating  turbulence  are  buoyancy  and  wind  shear.  Monin  and  Obukhov 
(1954)  used  this  concept,  combined  with  dimensional  analysis,  to  derive  the  relationship 
between  the  mean  wind,  temperature,  and  specific  humidity  profiles  to  surface  turbulent 
fluxes.  The  resultant  non-dimensional  gradients  can  be  expressed  as: 

Kz  du  ,  .  z  , 

— = 

Ut  oz  L 


Kz  do  ,  .  z  . 


(5) 


q*  oz  L 

where  von  Karman’s  constant /f  «0. 35-0. 40  derived  from  measurements,  L  is  the  Monin- 


Obukhov  length:  L  =  — ^  ,  g  is  gravity,  z  is  height,  and  <!>„(—)  and  (t\^(—)  are  non- 

Kgw’ O'  L  L 

dimensional  universal  functions  that  account  for  the  thermal  stability  of  the  surface  layer. 


The  left  hand  side  of  the  equations  is  the  non-dimensional  vertical  gradient  of 
mean  wind,  mean  potential  temperature,  and  mean  specific  humidity.  The  right  hand  side 
of  the  equations  indicates  that  these  mean  gradients  are  functions  of  height  non- 
dimensionalized  by  L,  denoting  the  effects  of  thermal  stability.  When  the  surface  layer  is 
neutral,  the  universal  function  for  momentum  is  equal  to  one.  The  other  universal 
function  values  have  been  determined  by  experiment.  Businger  et  al.  (1971)  derived  their 
widely  used  land-based  universal  functions  in  the  famous  1968  Kansas  experiment. 
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The  MOST  non-dimensional  functions  provide  the  link  between  bulk,  mean 
profiles  and  turbulence  measurements.  This  relationship  developed  by  Monin  and 
Obukhov  serves  two  primary  functions.  First,  MOST  allows  researchers  or  models  to 
calculate  surface  layer  turbulent  fluxes  using  of  mean  profile  measurements.  This 
procedure  is  especially  important  when  fast  turbulence  direct  sampling  is  not  possible. 
MOST  also  allows  researchers  to  calculate  mean  wind  at  a  level  at  which  they  did  not 
take  direct  observations.  For  example,  many  flux  parameterization  equations  require 
mean  winds  at  10m.  By  using  the  non-dimensional  gradient,  known  mean  values  at  one 
level,  and  surface  flux  values,  researchers  can  compute  the  wind  at  10m.  Finally,  if 
observations  are  available  at  only  one  level  and  flux  is  unknown,  both  the  gradient  and 
the  surface  turbulent  fluxes  can  be  derived  by  using  the  surface  roughness  height  and 
assumed  mean  properties  at  the  roughness  length.  The  surface  roughness  height  ( )  is 
the  height  at  which  mean  wind  is  zero  (Charnock  1955).  Similarly,  the  air  temperature  at 
the  surface  roughness  height  is  set  equal  to  SST,  and  the  specific  humidity  over  salt  water 
at  the  surface  roughness  height  is  defined  as  98  percent. 

MOST  does,  however,  have  limitations  and  researchers  continue  to  refine  the 
universal  function  values  and  debate  the  conditions  under  which  MOST  is  valid  (Token 
2006).  MOST  assumes  horizontally  homogenous  conditions  that  may  not  hold  true  near 
the  coast  or  near  a  front  (Cheney  2011).  Additionally,  MOST  assumes  that  the  constant 
flux  layer  is  in  equilibrium  with  the  surface  which  may  not  hold  true  over  rough  surfaces 
(Token  2006).  For  naval  and  marine  applications,  there  are  also  questions  as  to  whether 
the  Businger  universal  functions  derived  over-land  are  also  valid  over  water. 

D,  SURFACE  FLUX  PARAMETERIZATION 

Surface  flux  algorithm  calculates  surface  turbulent  fluxes  using  bulk,  mean 
measurements  because  bulk,  mean  quantities  are  more  readily  sampled  or  forecast  by 
numerical  models.  Per  Fairall  et  al.  (2003),  turbulent  flux  parameterizations  usually  use  a 
bulk  transfer  coefficient.  For  example,  the  flux  of  a  variable  x  would  be 
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where  is  the  bulk  transfer  coefficient,  S  is  the  mean  wind  speed  relative  to  the  ocean 

surface,  and  AX  is  the  difference  in  the  mean  value  of  x  from  the  sea  to  the  air.  The 
standard  bulk  expressions  for  the  primary  air-sea  surface  fluxes  are  then 

T  =  pC,S(U-U) 

H^=pCf,S(T^-e)  (7) 

Hi  =pL^C^S(q-q) 

where  Q ,  Q ,  and  are  the  total  transfer  coefficients  for  wind  stress,  sensible  heat,  and 

latent  heat,  respectively;  the  subscript  ‘s’  represents  the  value  at  the  air-sea  interface.  The 
representation  of  the  drag  and  exchange  coefficients  are  the  key  of  this  parameterization 
algorithm. 

Building  on  the  foundation  of  MOST  and  the  Businger  universal  function 
parameterizations,  the  Tropical  Ocean  Global  Atmosphere  Coupled-Ocean  Atmosphere 
Response  Experiment  (TOGA  COARE)  added  ocean  based  observations  to  the  flux 
database  and  refined  the  bulk  flux  parameterization  (Eairall  et  al.  1996b).  The  researchers 
used  these  simultaneously  collected  bulk  and  flux  measurements  to  refine  bulk 
algorithms  and  the  associated  transfer  coefficients,  roughness  parameter,  and  universal 
functions.  Their  modifications  attempted  to  account  for  light  winds,  cooling  effects  of 
precipitation,  and  convection  (Eairall  et  al.  1996b).  The  COARE  parameterizations  and 
their  subsequent  revisions  are  currently  widely  used  by  air-sea  researchers  and  coupled 
forecast  models  (Eairall  et  al.  2003).  The  analysis  presented  in  this  thesis  will  utilize  the 
COARE  algorithm  to  obtain  surface  fluxes  and  near-surface  mean  profiles. 

E.  UNCERTAINTIES  IN  SURFACE  FLUX  PARAMETERIZATION 

According  to  Bao  et  al.  (2002),  there  are  two  main  sources  of  uncertainty  for  sea- 
air  fluxes.  The  first  type  of  uncertainty  derives  from  uncertainty  in  the  observational  data 
and  measurements,  themselves.  The  second  major  source  of  uncertainty  results  from 
questions  of  when  MOST  is  valid.  In  particular,  there  remain  questions  about  how  well 
MOST  performs  in  conditions  with  very  high  or  low  winds,  waves  and  swells,  and 
different  surface  layer  stabilities. 
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1. 


Instrumentation  or  Data  Errors 


Measurements  of  the  bulk  and  turbulent  near  surfaee  quantities  eontinue  to  pose 
ehallenges  in  terms  of  how  the  observations  are  gathered,  used,  and  analyzed.  Tsvang  et 
al.  (1973)  note  the  importanee  of  instrument  intereomparisons  to  ensure  near  surfaee 
measurement  aceuraey.  In  their  series  of  turbulenee  measurement  experiments,  Tsvang  et 
al.  (1973)  note  that  “the  influenee  of  instrumentation  and  analysis  on  the  results  is  often 
subtle  and  complieated  so  that  it  may  not  be  readily  detected,  let  alone  removed.”  For 
example,  different  research  groups  have  obtained  different  values  for  universal  functions 
and  parameterizations.  These  differences,  however,  may  result  from  differences  in 
measurement  method  or  quality,  rather  than  inadequacies  of  the  underlying  theory 
(Tsvang  et  al.  1973).  Instrument  intereomparisons  help  to  eliminate  some  of  these 
uncertainties  and  potential  sources  of  error. 

SST  measurements  are  very  difficult,  if  not  impossible,  to  obtain  exactly  at  the 
interfacial  surface.  In  addition  to  flow  disturbance  caused  by  the  instruments  and 
observation  platforms,  the  warm  layer  and  cool  skin  effect  described  in  Fairall  et  al. 
(1996a)  affects  the  ability  to  obtain  accurate  interfacial  temperature  measurements.  The 
warm  layer  results  when  solar  radiation  warms  the  upper  few  meters  more  than  it  heats 
the  deeper  mixed  layer.  The  cool  skin  results  from  a  combination  of  net  longwave 
radiation  and  heat  flux  that  cools  the  upper  few  millimeters  mare  than  it  cools  the  layer 
beneath  it.  Thus,  depending  on  where  in  the  water  column  SST  is  measured,  in  situ 
measurements  may  need  to  be  adjusted  for  further  analysis  or  use  in  parameterization 
equations. 

2,  High  or  Low  Wind  Conditions 

There  are  questions  as  to  whether  MOST  and  the  algorithms  derived  from  that 
theory  remain  valid  in  very  high  or  very  low  wind  conditions.  MOST  and  the  bulk 
aerodynamic  algorithms  derived  from  MOST  show  that  the  drag  coefficient,  Cd,  should 
increase  with  wind  speed  over  the  water  (Zhu  and  Furst  2013).  In  very  low  winds, 
however,  observations  show  that  Cd  actually  increases  with  a  decrease  in  wind  speed 
(Zhu  and  Furst  2013).  In  very  high  winds,  spray  conditions  also  affect  heat  fluxes.  Total 
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sensible  and  latent  heat  fluxes  are  actually  composed  of  interfacial  fluxes,  which  are 
accounted  for  in  COARE,  and  fluxes  on  the  surface  of  sea  spray  droplets  which  are  not 
accounted  for  in  COARE  (Andreas  et  al.  2008).  Eor  wind  speeds  over  10  ms'\  the  spray 
fluxes  become  a  significant  fraction  of  the  overall  flux  and  cannot  be  neglected  (Andreas 
et  al.  2008). 


3.  Wave  Boundary  Layer  and  Effect  of  Waves 

According  to  Edson  and  Eairall  (1998),  “the  use  of  overland  measurements  to 
infer  surface  fluxes  over  the  open  ocean  raises  questions  about  the  universality  of  these 
relationships.”  Specifically,  Edson  and  Eairall  (1998)  theorize  that  wave-induced  effects 
in  the  wave  boundary  layer  (WBL)  can  create  conditions  where  turbulent  fluxes  do  not 
dominate  total  fluxes.  Rutgersson  et  al.  (2001)  explain  that  the  WBE  is  the  layer  at  the 
ocean  surface  that  is  directly  affected  by  waves,  and  in  which  the  surface  roughness 

height,  ,  changes  as  the  waves  evolve.  Charnock  (1955)  derived  a  relationship  for  the 
surface  roughness  length  as 


where  is  a  constant  called  the  Charnock  parameter;  m,  is  frictional  velocity,  and  g  is 

the  gravitational  acceleration.  As  winds  increase,  waves  and  the  surface  momentum  flux 
in  the  WBL  also  increase.  Thus,  per  the  Charnock  (1955)  equation,  the  surface  roughness 
height  in  the  WBL  is  dependent  on  the  wind  and  wave  fields.  According  to  Rutgersson  et 
al.,  “As  the  waves  evolve,  the  surface  roughness  changes,  so  the  flow  over  the  surface  is 
coupled  to  the  evolution  itself’  (2001).  The  so-called  sea-state  dependent 
parameterizations  attempts  to  account  for  this  changing  surface  roughness  and 
relationship  to  wave  age  and  wave  steepness.  Edson  and  Eairall  (1998)  do  conclude, 
however,  that  MOST  is  valid  over  the  ocean  for  measurements  and  calculations  above  the 
WBL. 
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F. 


SURFACE  LAYER  EM  PROPAGATION  AND  DUCTING 


Another  important  application  of  quality,  near  surface  mean  measurements  is  for 
use  in  EM  propagation  prediction.  EM  propagation  paths  depend  on  the  refractive 
structure  of  the  surface  layer  as  determined  by  vertical  temperature  and  humidity 
gradients.  Modified  refractivity,  M,  describes  this  structure  while  also  taking  into  account 
the  curvature  of  the  earth.  The  modified  refraction  index,  M,  can  be  given  by: 


M  ^N  +  0.157z^ 


4810e 

T 


)  +  0.157z 


(9) 


where  N  - 


4810e 

T 


is  the  radio  refractivity  index;  p  is  the  total  atmospheric 


pressure;  T  is  the  air  temperature  in  Kelvin;  and  e  is  the  water  vapor  pressure  in  millibar. 
There  are  four  general  propagation  path  types  defined  by  the  vertical  gradient  of  M:  sub¬ 
refraction,  standard  refraction,  super-refractions,  and  ducting.  Eigure  2  shows  the  four 
propagation  types  and  their  respective  modified  refractive  gradients.  The  ducting 
conditions  in  which  rays  are  bent  downward  or  “trapped”  when  dM/dz  <0,  are  of 
particular  interest  for  Naval  applications.  When  EM  waves  become  trapped  in  the  ducting 
layer,  the  duct  acts  as  a  waveguide  and  leads  to  decreased  propagation  loss  and  extended 
radar  ranges. 
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dA//dz  =  0 
^-'dW/dz  - 157 
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X 

X 


Sub-relraclion 


Figure  2.  Categories  of  radar  propagation  based  on  refractive  index  (from 

Turton  et  al.  1988). 


Evaporation  ducts  are  a  specific  type  of  surface  ducts  that  form  over  the  ocean 
due  to  a  steep  decrease  in  humidity  or  increase  in  temperature  with  height  (Ziemba  2013). 
Evaporation  ducts  often  occur  in  areas  where  there  is  a  temperature  inversion  combined 
with  a  significant  decrease  in  moisture  as  height  increases  (Ziemba  2013).  Eigure  3(f) 
shows  the  Mprofde  resulting  in  a  typical  evaporation  duct. 
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EXAMPLES  OF  MODIFIED  REFRACTIVITY  PROFILES 


MODIFIED  REFRACTIVITY  M 

Figure  3.  Typical  M  profiles  resulting  in  ducts  (from  Babin  et  al.  1997). 

The  strength  of  the  evaporative  duct  is  related  to  the  evaporation  duct  height 
(EDH).  The  EDH  is  defined  as  the  height  where  the  vertical  gradient  of  the  modified 
refractivity  changes  from  negative  to  positive.  As  Paulus  (1989)  explains,  the  EDH  “is 
not  a  height  below  which  an  antenna  must  be  located  to  give  extended  ranges  but  rather  a 
parameter  indicative  of  the  strength  of  the  evaporation  duct.”  Typical  ranges  for  the  EDH 
are  0-40m.  A  high  EDH  indicates  a  strong  evaporative  duct  that  can  trap  lower 
frequencies,  while  a  lower  duct  indicates  a  weak  evaporative  duct  that  can  only  affect 
higher  frequencies  (Paulus  1989). 

Although  ideally  one  would  measure  refractivity  and  EDH  directly,  it  is  not 
feasible  on  a  large  scale.  Therefore,  researchers  have  developed  models  based  on  bulk 
measurements  and  MOST  to  calculate  the  EDH  (Babin  et  al.  1997).  The  primary  factors 
that  determine  the  EDH  are  SST,  air  temperature,  relative  humidity,  wind  speed,  and  the 
height  of  the  measurements.  The  EDH  models  are  sensitive  to  the  bulk  measurement 
inputs,  especially  in  near-neutral  stability  conditions  (Edson  et  al.  1999).  Eigure  4  shows 
the  output  of  the  NPS  EDH  model  developed  by  Erederickson  in  2007.  In  near-neutral 
conditions  (air  temperature  approximately  equal  to  the  underlying  SST),  small 
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differences  in  wind  speed  led  to  veiy  large  changes  in  EDH.  Table  1  lists  the  Paulus 
(1989)  bulk  measuiement  precision  requuements  for  EDH  calculations. 


Figure  4.  Near'  neutral  condition  sensitivity  to  mean  measmement  inputs  (from 

Twigg  2007). 


Parameter 

Type 

Precision 

Wind  speed 

input 

1  kt 

Air  Temp 

input 

0.1  C 

RH 

input 

1% 

SST 

input 

0.1  C 

EDH 

output 

0.1  m 

Table  1 .  Bulk  measmement  precision  requirements  for  EDH  calculations 

(after  Paulus  1989). 
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G.  NEAR  SURFACE  MEASUREMENT  PLATFORMS 

While  platforms  such  as  ships,  fixed  buoys,  and  flux  buoys  already  do  exist  for 
collecting  near  surface  observations,  each  presents  challenges  or  limitations.  Measuring 
air-sea  wave  and  flux  directly  has  proven  difficult  for  many  practical  reasons.  Per  Edson 
et  al.  (1998),  “problems  usually  arise  from  three  sources:  platform  motion,  flow 
distortion,  and  environmental  factors  unique  to  the  ocean.”  Ship-based  measurements  are 
usually  obtained  from  an  elevated  level  above  the  sea  surface,  and  further  suffer  from 
flow  distortion  introduced  by  ship  superstructures.  As  described  above,  satellite  remote 
sensing  introduces  too  many  assumptions  imbedded  in  data  retrieval.  NOAA  NDBC 
buoys  provide  only  mean  measurements,  are  costly  to  maintain,  and  are  fixed  in  location. 
The  most  ideal  platform  for  obtaining  flux  measurements  are  flux  buoys,  which  are 
normally  large  and  heavy  (-2000  lb),  however.  Deployment  and  recovery  of  these  buoys 
is  difficult  and  costly.  Additionally,  these  buoys  usually  require  complex  mooring  in  deep 
waters. 

Taking  advantage  of  recent  development  in  sensor  miniaturization,  especially  in 
the  inertial  motion  sensing,  the  Naval  Postgraduate  School  (NPS)  Meteorology 
Department  developed  the  Marine- Air- Sea-Flux  (MASFlux)  buoy  (Zuniga  2013).  The 
MASFlux  buoy  is  an  instrumented  spar-buoy,  which  is  easy  to  deploy  and  recover,  and 
does  not  distort  the  flow  or  temperature  characteristics  of  the  surrounding  water  (Figure 
5).  This  system  measures  turbulent  perturbations,  mean  wind  and  thermodynamic 
variables  at  multiple  levels,  2-dimensional  surface  wave  spectra,  and  upper  ocean 
temperature  simultaneously.  The  small  buoy  is  not  suited  for  deployment  longer  than  a 
few  days,  however,  and  may  not  be  operable  in  harsh  conditions  such  as  strong  winds  (> 
15  ms'^)  and  high  seas.  A  detailed  description  of  the  NPS  MASFlux  system  can  be  found 
in  the  thesis  work  by  Zuniga  (2013). 
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Figures.  MASFlux  Buoy  deployment. 

In  addition  to  the  limitations  of  each  of  the  systems  described  above,  the  current 
instrument  suites  do  not  provide  truly  global  coverage  for  near  surface  observations.  The 
World  Meteorological  Organization’s  (WMO)  Global  Observing  System  (GOS)  relies  on 
approximately  11,000  stations  to  provide  meteorological  measurements  at  least  every 
three  hours.  Many  of  these  stations,  however,  and  land  based  and  there  remain  many  data 
sparse  regions,  especially  in  the  open  ocean  and  in  the  Southern  Hemisphere.  For  forecast 
models  in  particular,  errors  that  arise  in  data-sparse  regions  are  difficult  to  identify  or 
diagnose.  Not  only  do  these  coverage  gaps  reduce  accuracy  in  those  regions,  but  errors 
may  propagate  downstream  to  data-rich  regions  as  well  (Mathur  et  al.  1993). 
Additionally,  climate  change  researchers  require  accurate  ocean-atmosphere  flux  data  and 
there  are  only  limited  observations  in  polar  regions  where  heat  and  momentum  exchanges 
are  known  to  be  particularly  high  (WHOI  2011). 
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H.  WAVE  GLIDER  /  SHARC  OVERVIEW 


Liquid  Robotics,  Inc  (LRI),  a  Sunnyvale,  California  (CA)  based  eompany, 
introdueed  the  Wave  Glider  in  2007.  The  Wave  Glider  is  an  autonomous,  wave-propelled 
unmanned  surface  vehiele  (USV)  with  a  two-body  design.  When  used  by  the  Navy,  the 
Wave  Glider  is  known  as  the  Sensor  Hosting  Autonomous  Remote  Craft  (SHARC).  The 
SHARC  derives  its  power  from  wave  and  solar  energy  and,  therefore,  it  provides  a  long- 
endurance  and  long-range  platform  for  a  variety  of  applications. 

The  SHARC  consists  of  two  main  portions,  an  upper  “float,”  whieh  rides  on  the 
surface  of  the  water,  and  a  lower  “sub,”  whieh  is  tethered  approximately  seven  meters 
below  the  “float”  via  a  flexible  umbilical  (Figure  6).  The  float  and  sub  work  in  tandem 
eonverting  the  wave  aetion  energy  in  the  ocean  into  forward  propulsion.  As  the  float  rises 
and  falls  with  the  wave  motion  at  the  surface,  the  sub  portion  which  resides  in  relatively 
stationary  water  seven  meters  below  the  surface,  uses  rotating  “wings”  to  eonvert  some  of 
the  vertical  motion  into  forward  thrust  (Figure  7).  In  one  to  three  foot  seas,  the  SHARC 
maintains  an  average  speed  of  1.5  kts  and  up  to  one  year  of  endurance  (Hine  et  al.  2009). 


Figure  6. 


Illustration  of  main  eomponents  of  SHARC  (from  Roehholz  2010). 
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Reduced  wave  motion  at  depth 

Figure  7.  Illustration  of  how  the  wave  glider  converts  wave-powered  vertical 
motion  into  forward  motion  (from  LRI  2012a). 

Also  enabling  the  SHARC’s  long  endurance,  is  its  use  of  solar  power  and 
rechargeable  batteries  to  power  its  communications,  navigation  system,  and  additional 
payloads.  On  average,  the  SHARC  can  provide  10  watts  of  continuous  power  to  its 
modular  payloads  (Hine  et  al.  2009).  Communicating  with  the  SHARC  via  an  Iridium 
Satellite  Modem  and  shore-based  server,  an  operator  can  monitor  the  status  of  the  wave 
glider  and  instruments,  receive  small  packets  of  data,  and  input  course  adjustments.  The 
operator  can  also  prioritize  power  usage,  as  necessary,  by  dictating  which  sensors  are  on 
or  off.  All  of  these  aspects  have  proven  important  when  the  SHARC  is  operating  in 
hazardous  weather  conditions  (Griffith  et  al.  2012). 

Originally  intended  as  a  platform  for  listening  to  whales  (Rochholz  2012),  the 
Wave  Glider  can  also  be  used  as  a  sea-going  platform  to  make  measurements  of  the  air- 
ocean  environment.  In  fact,  one  of  the  most  attractive  features  of  the  SHARC  is  its 
capacity  to  deploy  different  payloads  based  on  the  mission  objectives.  The  SHARC  float 
contains  modular  electrical  interfaces  and  two  open  bays  (Figure  8)  for  additional 
payload  sensors  and  batteries. 
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Figure  8.  Illustration  of  the  Wave  Glider’s  modular  payload  dry  boxes  (from 

LRI  2012a). 

In  addition  to  the  open  payload  bays,  LRI  offers  both  the  basie  “METOC”  and  the 
“METOC  Plus”  models  of  the  SHARC  with  built-in,  integrated  weather  sensors  (Table 
2).  The  basic  LRI  METOC  Wave  Glider  model  uses  the  AIRMAR  PB  200  to  gather  wind 
speed  and  direction,  air  temperature,  and  barometric  pressure.  The  expanded  ERI  model 
(METOC  Plus)  adds  wave  and  water  sensors.  The  Datawell  MOSE-G  GPS-based 
directional  wave  sensor  allows  the  SHARC  to  measure  2-dimensional  wave  spectra,  wave 
height,  propagation  direction,  and  period,  while  the  Sea-Bird  GPCTD  allows  the  SHARC 
to  measure  water  temperature  and  salinity.  Most  ERI  gliders  acquired  by  Navy 
laboratories  and  centers  (e.g.,  NAVO  and  SPAWAR)  for  meteorology  and  oceanography 
applications  are  the  METOC  Plus  model.  Although  flux  instrumentation  is  not  part  of  the 
METOC  or  METOC  Plus  package,  the  SHARC  also  has  great  potential  as  a  platform  for 
collecting  surface  flux  measurements.  The  customizable  payload  bays  of  the  SHARC 
mean  that  collecting  METOC  data  and  performing  other  missions  are  by  no  mean 
exclusive  tasks. 
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Sensor 

Model 

Description 

METOC 

METOC 

Plus 

Airmar  FB200WX 

Weather  information  (air 
temperature,  wind  speed  & 
direction,  and  barometric  pressure) 

X 

X 

Weather  Station 

© 

Teledvne  RDI 

Measurement  of  ocean  currents  at 
the  surface  and  in  bins  to  the 
sensor’s  maximum  depth 

X 

X 

Acoustic  Doppler 

Current  Frofiler 

lADCF) 

Datawell  MOSE-G 

Wave  height,  direction,  period  and 
spectrum 

X 

GFS-based 

Directional  Wave 

Sensor 

@ 

Sea-Bird  GFCTD 

Water  salinity  analysis  through 
collection  of  conductivity, 
temperature  and  density;  dissolved 
oxygen  measurement 

X 

+  DO 

Table  2.  Table  Sensor  and  Instruments  for  LRI  METOC  gliders 

(after  LRI  2012b). 


1.  Airmar  PB200WX 

The  NFS  SHARC,  “Mako,”  is  a  basic  LRI  METOC  model  and  arrived  at  NFS 
equipped  with  the  Airmar  FB200WX  Weather  Station  for  near  surface  meteorological 
measurements.  Figure  9  shows  the  Airmar  FB200WX  mounted  on  a  1.12  m  mast  at  the 
center  of  the  Mako  float  (Figure  9).  A  blowup  of  the  Airmar  FB200WX  is  shown  in  the 
same  figure. 
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Figure  9.  Airmar  PB200WX  weather  station  mounted  on  the  SHARC. 


The  PB200WX  weather  station  measures  air  temperature,  apparent  wind  speed 
and  direetion,  and  barometrie  pressure.  Using  an  internal  eompass  and  GPS,  the  PB200 
ean  also  ealeulate  true  wind  speed  and  direetion  (PB200  teeh  manual).  Table  3  shows  the 
manufaeturer  speeifieations  of  the  Airmar  PB200. 

LRI  preeonfigures  the  Airmar  for  direet  integration  into  the  SHARC  eentral  data 
management  system.  Per  this  eonfiguration,  the  Airmar  PB200  to  sends  data  at  speeified 
intervals  to  the  SHARC.  At  a  one  seeond  interval,  the  PB200  sends  the  wind,  pressure, 
temperature,  and  vehiele  heading  messages.  At  a  10  seeond  interval,  the  PB200  transmits 
its  GPS  position  and  GPS  time.  The  SHARC  then  averages  the  data  over  10  minute 
intervals  and,  transmits  the  averaged  observation  to  the  WGMS  server.  The  raw  Airmar 
data,  however,  are  neither  stored  nor  transmitted  to  the  WGMS  server  (Kirklin  2014, 
personal  eommunication). 
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SPECIFICATIO\S 


Wnd  O  knatt  to  80  knatt  p  MPH  to  92  MPH7 

Wnd  Awrfi/t/ort.  OJ  Mors  pi 

lAIntf  Spmd  Aoattocv  9  OT  t>  S5T  (32T  t>  12271  no  fimb^ration*: 
Low  k¥nd  Speed: :  0  knots  to  10  knots  {0  MPH  to  11.2  MPH}: 

WAS  effofof  1  ik/ior  fS.l  MPHj  ♦iOfc  of  fHidino 
Hiqh  \^nd  Speeds :  10  knots  to  SO  knots  fll.2  MPH  to  S2  MPH}: 

WAS  ttrot  092  Anort  Q.S  MPHJ  ot  S9b  WAS. 
ivhklwt4‘/  ts  gr^ef 

Wnd  Spmd  Mamsc  y  In  vm  oonMom  *  *:  5  knots  (S.  7  MPH)  WAS 
Wnd  Ometon  Aoncr#.  cr  to  iecr 
tAtnd  Ometon  tkmoiuton.  QV 

MM  OmSon  Jtoaswey  #  OT  »  SSX  02T  »  132^  no  pmSpitoSan*. 
Low  Wind  Speeds.  Akms  tolOknts  f4  6MPHtoll.Sh$H  S'  WAStwksd 
High  Wnd  Speeds  :  »10  knots  (>11.2  MPH}~2*  .PMS  t,'dcc: 

MM  Ometion  Aacuroey  In  kwr  oonti$ons^: 

>8  knots /92  MPH)  S'  WAS  tyfliai) 

Carpass  Acoitocy: 

1*  sror>f  kwotStp  oomfocy 
2'  dynorrU  kwadnq  occutocy 
fkmof  Totn:  O*  ro  ^  (»f  second 
fkdoof  TUtn  Aoaeocyi  V  pet  secoyyd 

fkmof  Turn  am  ouput  ipdate  fkm 
2He-Nh(£A  0(S3  (Ahntatle  upWWHd 
At^abie  up  to  20  Nr-AIMC4  2000 
AiA  end  Roit  Range/Accufocy:  ±9Cr  /  <V 
At  "rntpetoPate  Range  2S'C  to  SS'C  flS'F  to  ISI'P) 

At  iwrperawie  Rmeduton:  QVC  pVP) 

At  Rnpmteute  Aixxeatp.  H'C^LS'Pf  9>4  A/urs  ^46  MTH)  wind 
Rm>nnettlc  Rmuun  Range: 

SSO  rmof  ro  iiso  rmot  os  inhg  to  34  inHg  SSO  hPo  to  SS.SO  Wo) 
aam>tnee/U  nmiiite  RauA/Uon:  Qlnttat  p£B9  inH^  QlWaf 
Bammettlc  Amst/te  Aoa/roey. 

±2  mPof  f±<X  OSS  i2  AAU  when  cottection  H  o^oHobie 

(PS  Aoc^Van  tkamey  3  m  (1(7)  lV4AS/raW5  fSSk.  of  the  time,  SA  off) 

Coating  'mrperotuee  Range  to  SS'C  (I3'f  to  I32'Fy 

Sigsply  Mbtwge.  9  VDC  to  16  WX 

Rgply  CUeene  <220  nA  AMght  2StS 

grane,  fQ7  Uh 

SmKOe  Baud  Rede  OMEA  02JB  bmefoee  Odd 
A80O  bps  icon  be  increased  to  38,400  bps  with  ocotmiond) 


Table  3.  Airmar  PB200WX  specifications  (after  Airmar  2011). 


2,  Existing  Datasets  of  SHARC/Airmar  Measurements 

Various  Wave  Glider  /  SHARC  user  groups  have  made  measurements  in  an 
attempt  to  evaluate  the  SHARC  Airmar  PB200WX.  This  research  is  especially  important 
because  the  Airmar  is  the  standard  weather  sensor  outfitted  on  the  SHARC.  Therefore, 
there  exists  a  potentially  vast  network  of  near  surface  measurements  obtained  by 
SHARCs  operating  over  the  world’s  oceans.  In  particular,  LRI  examined  their  Pacific 
Crossing  (PacX)  data  in  comparison  to  nearby  NDBC  buoys.  Naval  Oceanographic 
Office  (NAVO)  deployed  a  SHARC  in  the  vicinity  of  an  NDBC  buoy,  and  NPS  co¬ 
deployed  the  SHARC  and  the  MASFlux  buoy  for  weather  data  collection.  These  datasets 
are  to  be  analyzed  in  this  thesis. 
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a.  LRI  PacX  Measurements 

In  November  of  201 1,  LRI  deployed  four  Wave  Gliders  from  San  Francisco,  with 
two  heading  to  Australia  and  two  heading  for  Japan.  This  Pacific  crossing  served  as  both 
a  test  and  a  demonstration  of  Wave  Glider  capabilities.  The  data  collected  during  this 
crossing  is  known  as  the  “PacX”  data. 

During  this  deployment,  LRI  programmed  the  wave  gliders  to  orbit  several 
NDBC  buoys  so  that  they  could  compare  the  Airmar  /  SHARC  weather  measurements  to 
those  reported  by  the  NDBC  buoys.  On  February  18,  2012,  the  PacX  Wave  Gliders 
circled  buoy  46012  in  Monterey  Bay.  Similarly,  from  February  25-27,  2012,  three  of  the 
PacX  Wave  Gliders  orbited  within  eight  nm  of  NDBC  buoy  51000,  located  200nm 
northeast  of  Maui.  Griffith  et  al.  (2012)  focused  on  the  latter  dataset  for  their  report  on 
the  Wave  Glider’s  capability  as  a  weather  data  collection  platform. 

b.  Datasets  by  NA  VO  SHAR C  Deployments 

NAVO  evaluated  both  the  Airmar  and  the  SHARC  METOC  Plus  Model’s  GPS 
wave  sensor,  the  Datawell  MOSE  Datawell  GIOOO.  Eor  the  Airmar  evaluation  portion, 
NAVO  had  one  SHARC  circle  NDBC  buoy  51002,  near  Hawaii,  from  April  7-17,  2012. 
The  SHARC  made  a  hexagon  pattern  with  a  12km  radius  around  the  buoy.  This  area  also 
encompassed  one  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System  (COAMPS) 
model  grid  point,  and  one  U.S.  Navy  coastal  ocean  model  (NCOM)  grid  point.  Figure  10 
shows  the  track  of  the  NAVO  SHARC  during  this  evaluation  period.  NAVO  also 
provided  data  taken  intermittently  in  Monterey  Bay  from  September  1,  2012-December 
5,  2012.  Although  not  collecting  data  continuously,  the  SHARCs  collected  observations 
in  the  vicinity  of  NDBC  buoy  46092  for  a  total  of  approximately  12  days. 
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Figure  10.  NAVO  SHARC  track  as  it  orbited  NDBC  buoy  51002  in  April  2012 

(from  Wang  2012). 
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III.  SHARC  INSTRUMENTATION  AND  TESTING 


A.  DEVELOPMENT  OF  CUSTOMIZED  METOC  PAYLOAD  FOR  SHARC 

The  required  aceuraey  of  near  surface  measurements  necessitates  research  grade 
sensors  to  be  used  on  a  proven  platform  such  as  the  Wave  Glider.  Additionally,  as 
described  in  Chapter  II,  EDH  calculations  depend  on  humidity  measurements,  which  is 
not  available  on  the  default  METOC  Plus  sensor  package.  The  NPS  team  therefore 
decided  to  develop  a  sensor  package  for  the  SHARC  with  sensors  of  known  quality 
measurements.  This  system  was  mounted  on  a  1.1  m  mast,  aft  of  the  SHARC ’s  rear 
payload  bay  and  measures  wind  speed  and  direction,  barometric  pressure,  air 
temperature,  relative  humidity,  and  platform  position,  velocity,  heading,  pitch,  roll,  and 
six-axis  acceleration  rates  (Eigure  11).  A  separate  temperature  probe  also  provides  sea 
surface  temperature  (SST).  Table  4  lists  each  of  the  installed  sensors  and  their  measured 
variables.  Together,  the  NPS  sensor  suite  is  called  the  “NPS  Met.” 
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Additional  Data  Acquisition 

batteries  (in  System  (in  Aft  Bay) 

Forward  Bay) 

Figure  1 1 .  NFS  customized  METOC  sensor  suite  mounted  on  the  SHARC. 
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Sensor 

Measured  variables 

Temperature  and  Humidity 
Probe  Rotronic  Model 

MP100H 

Temperature  and  Relative 
humidity 

Vaisala  Weather  Transmitter 
WXT520 

Wind  speed  and  direction,  Barometric 
pressure 

Temperature 

Relative  Humidity 

Campbell  Scientific 
Temperature  Probe  Model 
109SS 

Sea  water  temperature 

Garmin  GPS16-HVS 

Position,  Velocity,  and 

Magnetic  declination 

True  North  Revolution 
Technologies  GS  Gyro 
Stabilized  Electronic 

Compass 

Heading,  Pitch,  and 

Roll 

VectorNav  VN-100  Rugged 
Accelerometer 

Angular  rates,  linear  accelerations,  and 
Magnetic  field  components 

Table  4.  NFS  customized  METOC  sensor  suite  and  measured  variables. 


A  CR3000  data  logger  encased  in  a  Pelican  model  1400NF  case  was  also  installed 
in  the  aft  payload  bay  of  the  SHARC.  The  dry  box  was  also  wrapped  in  Styrofoam  to  fit 
it  snugly  into  the  SHARC  payload  bay  (Figure  12).  All  of  the  additional  instruments 
described  above  connected  to  this  data  logger  as  the  data  acquisition  system  for  the 
sensor  suite.  The  CR3000  provides  approximately  7.5  hours  of  data  collection. 
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Figure  12.  Open  aft  payload  bay,  with  dry  box  and  data  logger  shown. 

The  new  sensor  payload  on  the  SFIARC  serves  two  purposes  during  the  NFS 
SFIARC  deployments  in  2014.  First,  the  new  sensor  suite  mirrors  the  sensors  used  on  the 
NFS  MASFlux  buoy,  a  system  that  had  been  previously  tested  and  proven  (Zuniga  2013). 
Co-deployment  of  the  SHARC  and  the  MASFlux  buoy  facilitates  a  direct  comparison  of 
the  datasets  collected  by  the  two  platforms.  This  comparison  ensures  we  are  able  to 
accurately  measure  the  mean  near-surface  properties  of  the  atmosphere  from  the  SHARC. 
Secondly,  it  provides  co-located  measurements  for  thorough  evaluation  of  the  Airmar 
FB200WX  with  known  quality  sensors.  This  configuration  allows  for  a  more  robust 
evaluation  of  the  Airmar  FB200WX  in  comparison  with  data  from  a  nearby  buoy  in 
previous  study. 
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B. 


NFS  SHARC  SENSORS 


The  following  paragraphs  describe  the  uses  and  design  of  the  sensors  NFS 
installed  on  the  SHARC.  Table  5  lists  manufacturer  stated  accuracy  and  ranges  of  the 
METOC  sensors,  GPS,  accelerometer,  and  compass. 


NPS  Sensor  and  Instrument  Specifications 

Temperature  and  Humidity  Proberotronic  Model  MPIOOH 

Parameter 

Temperature 

Humidity 

Range 

-40  to  60(C) 

Oto  100  (%RH) 

Accuracy 

0.03(K)(Ptl0O) 

0.3  (K)  (HygroClip) 

1.5  (%RH) 

Vaisala  Weather  Transmitter  WTX520 

Parameter 

Pressure 

Wind  Speed 

Wind  Direction 

Range 

600  to  1100  (hPa) 

Oto  60  (ms'^) 

Oto  360  (deg) 

Accuracy 

0.5(hPa)0to  30  C 
l(hPa)  -52  too,  30 to  60  C 

3%  at  10  ms'^ 

3  (deg) 

Resolution 

0.1  (hPa) 

.l(msh 

1  (deg) 

Campbell  Scientific  Temperature  Probe  Model  109SS 

Parameter 

Temperature 

Range 

-40  to +70  (C) 

Accuracy 

.49  at  -20  to  +70  (  C) 

Response  time 

36  (msec) 

Dip  angle 

+/-  80  (deg) 

Tilt  range 

+/-40(deg) 

Garmin  GPS  16-HVS 

Parameter 

Position 

Velocity 

Range 

+/-40(deg) 

+/-40(deg) 

Accuracy 

3  m 

0.1  RMS  (steady  state) 

PPS  Rate 

+/-1(HS) 

+/-l(ps) 

Repeatability 

+/-  0.2  (deg) 

+/-0.2(deg) 

VectorNav  VN-100  Rugged  Accelerometer 

Parameter 

Attitude 

Heading 

Acceleration 

Magnetic  Field 

Range 

+/-  90  (deg) 

+/- 180  (deg) 

+/-8(g) 

+/-  2.5  (Gauss) 

Accuracy 

0.5  (deg) 

2.0  (deg) 

+/-  0.05  (deg) 

+/-  0.05  (deg) 

Resolution 

0.05  (deg) 

Linearity 

0.5%  (FS) 

0.1%  (FS) 

Table  5.  Manufacturer  stated  accuracy  and  ranges  of  the  NFS  customizer 
METOC  sensors,  GPS,  accelerometer,  and  compass  package. 

1,  Rotronic  Model  MPIOOH  Temperature/Humidity  Probe 

This  instrument  includes  two  sensors  to  measure  relative  humidity  and 
temperature.  The  HygroClip  S3  samples  relative  humidity  while  the  Ft  100  samples  air 
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temperature.  Aeeording  to  Zuniga  (2013)  the  instrument  is  less  prone  to  errors  typieally 
eaused  by  solar  radiation  and  preeipitation.  Data  were  eolleeted  at  IHz. 


2,  Vaisala  Weather  Transmitter  WXT520 

The  WXT520  measures  and  ealculates  two-dimensional  wind  speed  and  direetion, 
and  barometrie  pressure.  The  sensor  uses  three  transdueers  and  ultrasound  teehnology  to 
determine  horizontal  wind  speed  and  direetion.  Importantly,  this  open  sensor  design  does 
not  impede  wind  flow  at  large  angles  of  platform  pitch  and  roll  as  compared  to  the 
Airmar  sonic  system  (Lind  2014,  personal  communication).  Data  were  collected  at  IHz. 

3,  Campbell  Scientific  Temperature  Probe  109SS 

A  CS  109SS  was  installed  in  the  aft  payload  bay  and  extends  through  a  small  hole 
in  the  underside  of  the  SHARC  hull  to  measure  SST.  The  probe  extends  approximately  1 
inch  from  the  hull  into  the  water.  On  the  third  day  of  field  testing,  three  small  Onset 
recording  thermometers  were  affixed  to  the  hull  to  provide  additional  SST  data  for 
comparison. 

4,  Garmin  GPS  16-HVS 

This  GPS  tracks  platform  speed  and  direction,  largely  used  to  calculate  true 
winds.  The  GPS  can  track  up  to  12  satellites  and  receives  one  second  navigation  updates. 
Additionally,  the  16-HVS  model  is  waterproof,  which  makes  it  well-suited  for  maritime 
use. 


5.  True  North  Revolution  Technologies  GS  Compass 

The  GS  model  compass  is  specifically  designed  for  platforms  operating  in 
dynamic  environments  such  as  rough  seas.  The  GS  uses  a  3-axis  solid  state  magnetometer 
to  provide  accurate  heading,  pitch,  and  roll  data.  Zuniga  (2013)  reported  reduced 
uncertainty  in  MASFlux  heading  data  after  adding  the  GS  compass. 
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6,  VectorNav  VN-100  Rugged  Accelerometer 

The  VN-100  Rugged  is  an  inertial  measurement  unit  (IMU)  and  Attitude  Heading 
Reference  System  (AHRS).  The  VN-100  combines  three-axis  accelerometers,  three-axis 
gyros,  and  3 -axis  magnetic  sensors  to  provide  three-dimensional  platform  orientation 
data.  According  to  Zuniga  (2013),  there  is  uncertainty  in  the  yaw  measurements.  Data 
were  collected  at  a  rate  of  5Hz. 

C.  FIELD  TESTS  OF  NFS  SHARC  INSTRUMENTATION 

The  NFS  Meteorology  Department  sought  to  further  evaluate  the  Airmar  in 
expanded  conditions  and  in  a  controlled  environment.  We  conducted  static 
intercomparisons  and  at-sea  testing  in  Monterey  Bay.  More  details  are  available  in  the 
below  sections. 

1.  Static  Intercomparison 

The  NFS  Meteorology  Department  conducted  a  comparison  of  the  Airmar  and  the 
new  instrument  package  in  a  controlled  environment.  Both  the  Airmar  and  the  new  sensor 
package  were  removed  from  the  SHARC  and  mounted  on  a  frame  as  shown  in  Figure  13. 
The  team  collected  data  on  the  roof  of  a  six  story  building  from  September  18  to  October 
6,  2014.  The  two  masts  were  then  moved  to  an  open  area  on  Fort  Ord  in  Marina,  CA 
from  October  6-24,  2014  in  order  to  expose  the  sensors  to  higher  wind  conditions.  From 
November  26  through  December  6,  2014,  the  static  test  resumed  atop  the  six  story 
building  to  confirm  some  of  the  initial  data.  These  static  datasets  allow  for  direct 
comparison  of  the  IHz  data  from  the  new  NFS  sensors,  with  the  raw  IHz  data  from  the 
Airmar  (as  opposed  to  the  10-minute  averaged  data  integrated  into  the  LRl  data 
management  system).  Additionally,  the  static  test  removes  the  platform  motion  from  the 
problem. 
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Figure  13.  Static  testing  of  the  Airmar  PB200WX  and  the  NFS  sensor  package 
a)  atop  a  six  floor  building  and  b)  at  Fort  Ord,  Marina,  CA. 

2,  At-Sea  Testing 

NFS  conducted  three  field  tests  of  the  SHARC  in  Monterey  Bay  on  July  29,  July 
30,  and  August  8,  2014,  respectively.  Not  only  did  the  SFIARC  host  the  new  sensor 
package  as  described  above,  the  SFIARC  was  also  co-deployed  with  the  NFS  MASFlux 
buoy.  The  SHARC  received  frequent  steering  updates  to  keep  it  within  2km  of  the 
MASFlux  buoy  throughout  the  five  hours  of  data  collection  each  day. 

a.  MASFlux  Configuration 

The  MASFlux  buoy,  with  its  near  surface  measurements  and  small  size,  is  most 
suitable  for  evaluating  the  SHARC  meteorological  and  upper  ocean  measurements.  The 
MASFlux  buoy  features  well  calibrated  wind,  temperature,  and  humidity  sensors  at 
multiple  levels  within  the  lowest  4  m  of  the  surface,  concurrent  wave  and  upper  ocean 
temperature  measurements,  and  turbulent  flux  measurements  using  direct  eddy 
correlation  method,  which  is  the  ‘golden  standard’  flux  measurement  compared  to  bulk 
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methods.  As  described  above,  the  MASFlux  buoy  consists  of  the  same  sensors  as  the 
SHARC:  the  Rotronic  humidity  and  temperature  probe,  Vaisala  weather  transmitter 
WXT520,  Campbell  Scientific  temperature  probe  109SS,  Garmin  GOS16-HVS,  True 
North  GS  Electronic  Compass,  and  VectorNav  VN-100  Rugged  accelerometer.  With  its 
similar  configuration  and  known  quality,  the  MASFlux  data  can  be  used  to  examine  the 
feasibility  of  using  the  current  METOC  glider  measurements  in  quantifying  various 
properties  at  the  air-sea  interface  from  wave  glider  measurements,  and  for  further  use  in 
surface  flux  parameterizations  and  forecast  models. 

Additionally,  the  MASFlux  buoy  hosts  the  Ultrasonic  Anemometer  R  M  Young 
Model  81000VRE.  This  sensor  provides  critical  measurements  for  obtaining  turbulent 
flux  calculations.  In  particular,  this  instrument  measures  three-dimensional  wind  velocity, 
and  the  speed  of  sound  at  a  rate  of  20  Hz.  For  details  of  the  MASFlux  instrumentation, 
refer  to  Zuniga  (2013). 

Figure  14  shows  the  MASFlux  buoy  and  SHARC  tracks  from  each  day.  Different 
weather  conditions  were  encountered  on  the  three  days  of  field  testing  and  different 
techniques  for  co-deploying  the  SHARC  and  MASFlux  were  tested.  On  both  July  29  and 
30,  2014,  winds  and  seas  were  calm.  On  August  8,  2014,  winds  again  were  calm  but  there 
were  1.5m  swells  and  confused  waves.  Table  6  shows  the  weather  conditions  for  each 
day  as  reported  by  nearby  NDBC  buoys. 
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SHARC  and  BuoyTracks 


Figure  14.  SFIARC  and  MASFlux  tracks  during  at-sea  field  testing. 


NDBC  Buoy  46236 
(Monterey  Canyon  Outer) 

29-Jul-14 

30-Jul-14 

8- Aug- 14 

Wave  Height  (m) 

0.86 

1.02 

1.61 

Dominant  Wave  Period  (s) 

8.33 

8.33 

8.33 

Avg  Wave  Period  (s) 

6.49 

5.77 

6.64 

SST  ( C) 

17.8 

17.7 

17.4 

NDBC  Buoy  46092  (MBARI) 

Wind  Speed  (m/s) 

3 

5 

7 

Table  6.  One  hour  averages  as  reported  by  NDBC  buoys  at  1102  PDT  on 

each  deployment  day. 


The  sensors  performed  well  each  day  with  a  few  exceptions.  On  July  30,  RA^ 
Fulmar  remained  too  close  to  the  MASFlux  buoy  while  the  buoy  was  deployed.  Frequent 
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backing  down  of  the  Fulmar  created  ocean  mixing  and  possibly  flow  disturbance.  The 
situation  was  corrected  for  future  testing  days  by  having  the  research  vessel  maintain  at 
least  .5  km  separation  from  the  SHARC  and  MASFlux  buoy  during  data  collection. 
Additionally,  extra  temperature  probes  were  added  to  the  SHARC  for  the  August  8 
deployment  to  gather  more  data  for  comparison. 

b.  Deployment  and  Recovery  Procedures 

All  deployments  and  recoveries  were  conducted  from  the  stem  of  the  (RA^) 
Fulmar,  a  66 ’9”  catamaran.  The  Fulmar  features  a  20001b  A-frame  that  the  crew  used  for 
hoisting  the  SHARC  into  and  out  of  the  water.  These  evolutions  are  dynamic  and 
challenging,  and  the  NFS  team’s  proficiency  in  deploying,  recovering,  and  operating  the 
SHARC  increased  on  each  deployment  day. 

In  particular,  deployment  and  recovery  procedures  of  the  wave  glider  were 
improved  and  codified  over  the  three  underway  days.  For  SHARC  recovery,  the  NFS 
crew  must  communicate  with  the  R/V  Fulmar  captain  to  intercept  the  SHARC  by 
maneuvering  the  Fulmar  so  that  the  Fulmar’s  fantail  is  on  the  glider’s  starboard  side. 
Then,  it  is  critical  that  the  captain  place  the  Fulmar  in  neutral  to  avoid  backing  over  the 
SHARC.  Attaching  fore  and  aft  taglines  to  SHARC  before  deployment  (and  leaving  them 
attached  throughout  that  day’s  underway)  also  greatly  improves  both  the  deployment  and 
recovery  processes.  Not  only  were  the  tag  lines  used  to  stabilize  the  SHARC  while  being 
lifted  by  the  A-frame,  but  the  tag  lines  also  allowed  the  Fulmar  to  stay  farther  away  from 
the  SHARC  during  recovery.  Finally,  during  data  collection,  the  Fulmar  must  maintain  at 
least  1000m  separation  from  both  the  SHARC  and  the  MASFlux  buoy  to  prevent  the 
Fulmar  from  disturbing  the  flow  characteristics  and  environmental  conditions  being 
measured. 
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IV.  NPS  MET  SENSOR  PERFORMANCE  AND  AIRMAR 

EVALUATION 


The  NPS  team  evaluated  the  new  NPS  Met  package  and  the  standard  SHARC 
weather  sensor,  the  Airmar  PB200,  against  known,  high-quality  environmental 
measurements.  First,  the  team  compared  the  NPS  Met  data  to  that  collected  from  the 
MASFlux  buoy.  Not  only  did  this  step  allow  the  team  to  determine  the  quality  of  the  NPS 
Met  suite,  but  it  also  allowed  the  team  to  determine  whether  the  SHARC  platform  adds 
distortion  to  the  near  surface  measurements  collected  from  the  SHARC.  This  evaluation 
aims  to  determine  the  suitability  of  both  the  SHARC  platform  and  the  NPS  Met 
suitability  for  air-sea  data  collection.  Additionally,  after  validating  the  NPS  Met 
instruments,  they  can  provide  a  set  of  known,  quality  measurements  against  which  to 
compare  the  Airmar  data.  The  Airmar  data  collected  by  NPS  then  are  compared  to  the 
NPS  Met  data,  and  the  Airmar  measurements  provided  by  NAVO  and  LRI  are  compared 
to  NDBC  Buoy  46092  in  Monterey  Bay. 

A,  NPS  MET  PAYLOAD  EVALUATION 

As  described  in  Chapter  II,  the  NPS  Met  sensors  mirror  those  used  on  the 
MASFlux  buoy.  The  following  sections  show  the  comparison  between  the  NPS  Met  and 
MASFlux  wind,  specific  humidity,  pressure,  and  temperature  measurements  collected  at 
sea.  In  general,  a  separate  time  series  was  plotted  for  each  variable,  and  for  each  day  of  at 
sea  testing.  The  scatter  plot  uses  the  combined  data  from  all  three  days  of  underway 
testing.  The  red  line  on  the  scatter  plot  represents  a  one-to-one  relationship  between  the 
two  data  sets.  Perfectly  correlated  and  unbiased  data  would  fall  on  that  line.  Both  the  time 
series  and  scatter  plot  graphs  show  a  10  minute  moving  average  of  the  data.  This  moving 
average  permits  more  direct  comparison  of  the  mean  values  from  the  different  sensor 
packages  or  platforms  and  smooths  the  datasets. 

1,  Wind  Speed  and  Direction  Measurements 

Wind  speed  measurements  for  each  testing  day  show  wind  speed  generally 
increasing  with  height,  which  is  consistent  with  MOST  (Figure  15).  On  each  day  of 
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testing,  the  winds  gradually  inereased  over  time  but  remained  light.  On  the  SHARC,  the 
NFS  Met  sampled  wind  speed  at  a  height  of  1.1  m.  This  height  falls  between  the 
MASFlux  wind  speed  sensor  levels  one  and  two.  Therefore,  we  expeet  the  NFS  Met  wind 
speed  to  fall  between  the  MASFlux  levels  one  and  two  wind  speeds.  Graphieally,  this 
means  the  magenta  line  (NFS  Met)  should  be  in  between  the  green  and  red  lines 
(MASFlux  levels  one  and  two,  respeetively).  With  minor  exception,  this  trend  does  hold 
true  on  all  three  days. 

Figure  16  shows  the  strong  correlation  between  the  NFS  Met  wind  speed  data  and 
the  MASFlux  level  one  on  all  three  days.  The  small  positive  bias  of  the  NFS  Met  is  due 
to  the  fact  that  the  NFS  Met  wind  speed  sensor  is  mounted  higher  than  the  level  one 
MASFlux  buoy  wind  sensor.  Also,  although  the  MASFlux  buoy  and  the  SHARC  were 
within  1.5  km  of  each  other  during  all  at-sea  testing,  they  were  not  exactly  collocated. 
This  small  spatial  variability  also  accounts  for  small  differences  between  the  two 
platforms’  measurements. 
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Figiue  15.  Time  variation  of  the  MASFliix  and  NFS  Met  wind  speed 

measmement  on  a)  July  29,  2014;  b)  July  30,  2014;  c)  August  8,  2014. 
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At-Sea  Test:  MASFlux  vs  NPS  Met  Wind  Speed 


Figure  16.  Comparison  of  the  level  one  MASFlux  wind  speed  measurements 
and  the  NPS  Met  wind  speed  measurements  during  at  sea  testing. 

Comparisons  of  the  wind  direetion  during  at-sea  testing  are  shown  in  Figure  17. 
Winds  remained  Southwest  to  Northwest  on  all  three  days  and  should  not  vary  over  the  4 
m  height  of  the  NPS  Met  and  MASFlux  sensors.  Therefore,  the  wind  direetion 
measurements  should  be  identieal  from  all  of  the  MASFlux  wind  sensors.  In  general,  the 
NPS  Met  wind  direetion  measurements  follow  the  same  trend  as  the  MASFlux 
measurements  but  with  some  variation.  The  scatter  plot  shows  generally  good  correlation 
(Figure  18)  between  NPS  Met  and  MASFlux  level  one  wind  direction  measurements. 
Wind  direction  measurements  from  July  30,  2014  show  slightly  decreased  correlation 
compared  to  the  other  two  days  of  at-sea  testing.  Again,  spatial  variation  accounts  for 
some  differences  in  the  two  datasets.  Additionally,  removing  the  movement  of  the 
SHARC  platform  from  the  NPS  Met  observed  winds  was  more  challenging  when  the 
SHARC  made  a  lot  of  turns.  On  day  two  of  testing  (July  30,  2014),  the  NPS  team  steered 
the  SHARC  in  a  tighter  pattern  to  stay  within  0.5- 1.0  km  of  the  MASFlux  buoy.  This 
tight  pattern,  increased  SHARC  turning,  and  more  difficult  data  processing  may  account 
for  differences  in  the  wind  direction,  as  well. 
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Figure  1 7.  Same  as  Figaire  15,  except  for  wind  direction. 


At-Sea  Teat;  MASFlux  vs  NPS  Met  Wind  DM 


Figure  18.  Same  as  Figure  16,  except  for  wind  direction. 
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2,  Specific  Humidity 

Figure  19  shows  the  time  series  plots  of  the  at-sea  speeific  humidity  that  was 
caleulated  using  the  relative  humidity  measurements  from  the  NFS  Met  and  MASFlux 
sensors.  On  July  29  and  July  30,  2104,  the  level  one  MASFlux  RH  sensor  went 
underwater  during  deployment  and  was  inoperable  for  the  remainder  of  those  two  days. 
Therefore,  only  MASFlux  levels  two  through  four  are  shown  for  July  29  and  30,  2014. 
Due  to  surfaee  evaporation,  RH  and  speeifie  humidity  deerease  sharply  with  height  as 
seen  in  the  lowest  -'20  m  of  the  atmosphere.  The  NFS  Met  samples  relative  humidity  at  a 
height  of  0.6  m  and  should  have  speeific  humidity  values  greater  than  those  measured  by 
the  MASFlux  level  two,  and  less  than  those  values  measured  by  level  one  of  the 
MASFlux,  when  available.  This  trend  is  observed  on  all  three  days  of  testing.  Figure  20 
shows  very  strong  correlation  between  NFS  Met  and  the  level  two  MASFlux  specific 
humidity.  The  small  bias,  approximately  0.15  gKg'^,  is  likely  caused  by  the  vertical 
gradient  and  spatial  variability. 
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At-Sea  Test:  MASFlux  vs  NFS  Met  q 


Figiue  20.  Same  as  Figure  16,  except  for  specific  humidity. 


45 


3. 


Pressure 


Figiue  21  is  the  time  series  plots  of  the  at-sea  NPS  Met  and  MASFlux  pressiue 
measurements  and  Figure  22  shows  the  scatter  plot  for  all  three  days.  There  is  nearly 
perfect  conelatiou  between  the  two  platforms  and  their  respective  instnrmeuts.  The  NPS 
Met  pressitre  sensor  has  a  bias  of  -0.14  mb  which  is  withirr  the  accirracy  limits  of  the  two 
sensors.  Tlierefore,  the  pressiue  measiuements  from  the  NPS  Met  and  the  MASFlux  buoy 
are  statistically  identical. 
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Figiue  2 1 .  Same  as  Figme  1 5,  except  for  pressiue. 
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At-Sea  Test;  MASFlux  vs  NFS  Met  Pressure 


Figure  22.  Same  as  Figure  16,  except  for  pressure. 

4,  Temperature 

Figure  23  shows  the  SST  and  air  temperature  measurements  from  the  NFS  Met 
and  MASFlux  sensors.  On  each  of  the  three  at-sea  days,  SST  was  warmer  than  air 
temperature,  denoting  unstable  thermal  stratifications.  The  air  temperature  decreased  with 
height  and  SST  decreased  with  increasing  depth,  as  expected.  The  NFS  Met  SST,  in 
general,  does  not  show  as  much  variability  as  the  MASFlux  measured  SST.  On  July  30, 
2014,  the  R/V  Fulmar  stayed  within  approximately  .25  km  of  the  MASFlux  buoy  and 
periodically  backed  down  to  maintain  position.  This  action  caused  a  lot  of  mixing  in  the 
area  of  the  MASFlux  buoy,  brought  cooler  water  to  the  ocean  surface,  and  introduced 
additional  variability  in  the  sampled  SST.  On  August  8,  2014,  there  was  high  swell  so  a 
greater  distance  between  the  SFIARC  and  MASFlux  buoy  had  to  be  used.  The  sharp 
gradient  in  SST  visible  on  that  day  may  be  an  ocean  front  and  the  distance  between  the 
two  platforms  could  explain  why  the  two  platforms  encountered  the  front  at  different 
times. 


For  air  temperature,  the  NFS  Met  samples  the  air  at  0.6  m  which  is  between  the 
MASFlux  levels  one  and  two.  On  Jul  29  and  July  30,  2014,  the  MASFlux  level  one 
temperature  sensor  malfunctioned  for  the  same  reason  described  above  in  section  2. 
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Therefore,  only  MASFlux  levels  two  through  four  are  shown  for  those  days.  Since  air 
temperature  decreases  with  height  in  the  observed  unstable  condition,  the  NFS  Met  air 
temperature  should  be  greater  than  those  measured  by  the  MASFlux  level  two,  and  less 
than  those  values  measured  by  level  one  of  the  MASFlux,  when  available.  This  trend  is 
observed  on  all  three  days  of  testing  with  a  few  deviations  likely  caused  by  the  spatial 
variability  described  above.  Figure  24  shows  strong  correlation  and  an  absolute  error  of 
only  0.05  degrees  Celsius  between  the  NFS  Met  air  temperature  measurements  and  the 
MASFlux  level  two  air  temperature  measurements. 
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Figure  23.  Same  as  Figure  15,  except  for  temperatme. 
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At-Sea  Test:  MASFlux  vs  NPS  Met  Temperature 


Figure  24.  Same  as  Figure  16,  exeept  for  air  temperature. 

5,  Solar  Heating  Effect 

The  data  above  show  that  the  effeet  of  the  platform  on  air  temperature  is  minimal 
and  sampling  error  from  the  NPS  Met  is  likely  from  other  sourees.  Solar  heating  of  the 
instrument  housing  is  known  to  affeet  the  aeeuraey  of  all  naturally  ventilated  air 
temperature  sensors,  sueh  as  the  Rotronie  model  used  on  both  the  MASFlux  and  NPS 
Met.  As  stated  in  Chapter  11,  air-sea  fluxes  and  near  surfaee  gradient  ealeulations  require 
high  quality,  aeeurate  mean  measurements.  Therefore,  while  this  solar  heating  will  not  be 
apparent  in  the  NPS  Met  and  MASFlux  eomparison  testing,  it  is  still  important  to 
eonsider  this  effeet  when  attempting  to  determine  temperature  aeeurately. 

An  R.  M.  Young  radiation  shield  was  used  with  the  Rotronie  MPlOO  to  reduee 
solar  heating  of  the  sensor.  Even  with  the  housing  shield,  however,  the  radiation  heating 
effeet  may  still  be  sizeable  in  some  eonditions  sueh  as  low  winds  and  elear  skies.  Figure 
25  shows  the  relationship  between  solar  heating,  the  natural  ventilation  rate,  and  the  sun 
angle.  In  general,  one  would  expeet  to  see  a  signifieant  solar  heating  effeet  in  low  wind 
eonditions  and  when  the  sun  is  lower  in  sky.  Of  note,  the  Airmar  does  not  have  an 
effeetive  radiation  shield  and  is  therefore  likely  more  suseeptible  to  heating  from 
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incoming  irradiation  than  the  NPS  Met  instrument.  This  point  will  be  diseussed  in  a  later 
section. 


DATA  SUMMARIZED  FROM  *COMf>ARlSON  TESTING  OF  SELECTED  NATURALLY 
VENTHATEO  RADIATION  SHIELDS*  BY  GERALD  C.  OIU.  SEPTEMBER  19B3. 


Figure  25.  Relationship  of  the  solar  heating  of  the  multi-plate  radiation  shield 
used  on  the  Rotronic  sensor  to  the  ventilation  rate  and  sun  angle  (from 

R.  M.  Young,  2012.) 


In  order  to  identify  the  signature  of  solar  heating  in  the  temperature  data,  the  NPS 
Met  and  a  solar  irradiance  sensor  eolleeted  data  simultaneously  at  Fort  Ord.  Figure  26 
shows  the  data  collected  on  October  9  and  October  10,  2014.  Oetober  9,  2014  is  a  typical 
sunny  day  as  evident  by  the  high  peak  value  of  solar  irradiance  and  smooth  irradiance 
eurve.  Figure  27  shows  the  data  collected  on  Oetober  17  and  Oetober  18,  2014  and 
represents  a  typical  cloudy  day  as  evident  by  the  lower  peak  solar  irradiance  value  and  a 
“fuzzy”  irradianee  eurve.  On  the  sunny  day,  there  are  two  peaks  in  NPS  Met  air 
temperature  whieh  correspond  to  the  period  of  steep  changes  in  solar  irradianee.  It  is  also 
possible  that  some  direct  sunlight  was  reaching  the  NPS  temperature  sensor  at  lower  sun 
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angles.  On  the  cloudy  days,  in  contrast,  there  are  no  extra  peaks  in  the  NFS  Met 
temperature  data  due  to  the  light  being  more  diffuse. 


Wind  Speed  09-0ct-2014-10-0ct-2014 


Figure  26.  Typical  measurements  of  NFS  Met  temperature,  incoming  solar 
irradiance,  and  wind  speed  at  Fort  Ord  during  the  intercomparison 
testing  on  a  sunny  day. 


Wind  Speed  17-0ct-2014-18-0ct-2014 


Figure  27.  Same  as  Figure  26,  except  for  a  cloudy  day. 
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6. 


NPS  Met  Mean  Measurement  Conclusions 


Overall,  the  NPS  Met  and  MASFlux  mean  measurements  are  highly  correlated 
and  display  only  small  biases  due  to  spatial  and  temporal  variability.  Table  7  summarizes 
the  statistical  comparison  between  the  NPS  Met  and  MASFlux  buoy  data.  The  SFIARC 
platform  does  not  introduce  significant  flow  distortion  for  the  mean  measurements  in 
comparison  with  the  small  and  drifting  MASFlux  buoy.  Therefore,  the  NPS  team  feels 
confident  of  the  NPS  Met  payload  system  and  will  use  it  as  the  “truth”  for  evaluating 
Airmar  data. 


MASFlux  (Levels  Indicated)  vs  NPS  Instruments  At  Sea  (Moving  10  Minute  Averages) 

Press  (mb) 

Temp  L2  (  C  ) 

wind  Sp  LI  (ms’’) 

RH  (%) 

Mean  Error 

0.14 

-0.08 

-0.31 

0.52 

Error  Std  Dev 

0.01 

0.06 

0.29 

0.32 

Abs  Error 

0.14 

0.08 

0.33 

0.55 

Corr  Coef 

1.00 

0.99 

0.97 

0.99 

Table  7.  Statistical  comparison  of  MASFlux  and  NPS  Met  measurements 

during  at  sea  testings. 


B.  NPS  TESTING  AND  EVALUATION  OF  THE  AIRMAR 

The  NPS  team  conducted  both  at-sea,  and  land-based  static  testing  of  the  Airmar 
PB200  weather  station.  Analysis  of  these  data  and  their  comparison  with  the  NPS  Met 
data  collected  concurrently  is  detailed  below. 

1,  NPS  Airmar  at  Sea  Testing 

The  SFIARC  transmits  10  minute  averaged  Airmar  weather  data.  To  evaluate  the 
at-sea  data,  I  averaged  10  minutes  of  NPS  Met  data  for  direct  comparison  with  the 
Airmar  data  from  the  SHARC.  The  Airmar  and  NPS  Met  were  separated  by  only  0.5  m 
on  the  SHARC  which  virtually  eliminates  spatial  and  temporal  variability  between  the 
two  datasets.  Therefore,  if  both  instruments  are  providing  high  quality  measurements,  the 
data  should  be  statistically  identical. 

Figure  28  is  the  time  series  plots  of  the  averaged  Airmar  and  NPS  Met  pressure 
data  collected  during  the  three  days  of  at-sea  testing.  While  the  Airmar  captured  the 
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overall  pressiue  trend  each  day,  it  did  not  handle  the  details  well.  The  Ainnar  tended  to 
luiderestimate  the  pressure  and  introduced  a  large  amount  of  variability  not  present  in  the 
NFS  Met  data.  Figme  29  shows  the  low  coixelation  between  the  Aiimar  and  NFS  Met 
pressure. 


20140729 


20140730 


20140808 


•-NPS  Met  (1.1  m) 
— Airmar  (1.1  m) 


Figure  28.  Time  variation  of  Aiimar  and  NFS  Met  pressure  measiuements 
dining  at  sea  testing  on  a)  July  29,  2014;  b)  July  30,  2014;  c)  August  8, 

2014. 
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At-Sea  Test:  Airmar  vs  NPS  Met  Pressure 


Figure  29.  Comparison  of  the  Airmar  pressure  measiuemeuts  versus  the  NPS 
Met  pressiue  measiuemeuts  dming  at-sea  testing. 

The  Airmar  performed  better  with  measming  wind  speed  and  wind  direction  at 
sea,  as  shown  in  Figme  30.  The  data  are  not  statistically  identical  but  they  are  well- 
correlated  as  shown  in  Figiue  31. 
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Figirre  30.  Time  variation  of  Airmar  and  NPS  Met  wind  speed  measiuemeuts 
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dmiiig  at  sea  testing  on  a)  July  29,  2014;  b)  July  30,  2014;  and  c) 
August  8,  2014;  and  time  series  plots  of  wind  diiection  measiu  ements 
on  d)  July  29,  2014  e)  July  30,  2014  and  f)  August  8,  2014. 


Figiue  3 1 .  Same  as  Figiue  29,  except  for  a)  wind  speed;  b)  wind  diiection. 


2.  NFS  Airmar  Static  Testing 

The  NFS  team  made  static  intercompaiison  using  the  Aiimar  and  NFS  Met  in 
order  to  evaluate  the  Aiimar’s  raw  1  Hz  data.  Figure  32  shows  the  time  variation  of  the 
raw  data  for  each  of  the  measiued  variables.  Red  indicates  Aiimar  measiuements  while 
blue  indicates  NFS  Met  measurements.  Just  as  with  the  at-sea  data,  the  Aumar  captiues 
the  pressiue  and  temperatme  trends  but  not  the  details,  and  introduces  additional 
variability.  The  Aiimar  performs  well  paiiiculaiiy  with  wind  measiuements. 

Of  note,  the  Aiimar  overestimated  the  pressiue  dming  the  static  test,  whereas  it 
imderestimated  the  pressiue  diuing  the  at-sea  testing.  The  10  minute  averages  of  the  static 
test  data  from  the  SHARC  data  management  computer  also  remained  higlier  than  the 
NFS  Met  measiuements,  indicating  no  additional  calibration  of  the  Aiimar  original  data 
by  the  LRI  data  system.  The  bias  of  about  5  mb  in  pressiue  seems  to  be  consistent  with 
previous  studies  of  LRI  (Griffith  et  al.  2012)  and  Naval  Research  Laboratory  (Wang  and 
Allard,  2012)  although  the  latter  were  both  from  measiuements  at  sea.  The  negative  bias 
fr  om  NFS  at-sea  testing  is  thus  puzzling  and  requires  more  at-sea  testing  to  frilly  address 
and  quantity. 
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static  Test,  Round  2 


29.5  30  30.5 

Day  of  November  2014  (UTC) 


Figure  32.  Time  variation  of  Airmar  1  Hz  measurements  during  statie  testing 
where  blue  and  red  indicates  NFS  Met  and  Airmar  measurements, 

respectively. 


The  Airmar  temperature  is  consistently  higher  than  the  NFS  Met  temperature 
during  the  daytime.  The  largest  difference  at  the  mid-day  can  be  as  high  as  5°C.  At  night, 
the  Airmar  temperature  is  slightly  higher  with  a  much  smaller  discrepancy  compared  to 
the  daytime  overestimated  values  (Figure  32).  This  diurnal  variation  clearly  indicates  the 
effect  of  solar  radiation  on  the  Airmar  temperature  measurements.  This  is  expected  as  the 
Airmar  temperature  sensor  is  housed  in  a  ~3  mm  diameter  void  with  plastic  surroundings. 
This  housing  was  not  designed  for  optimal  natural  ventilation  as  the  Rotronic  housing. 

Figure  33  shows  the  comparison  of  the  Airmar  pressure,  temperature,  and  wind 
speed  with  the  respective  NFS  Met  measurements.  To  further  study  the  effect  of  solar 
heating  on  the  Airmar  temperature,  I  separated  the  data  into  day  and  night  segments.  Data 
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collected  between  1 400-0  lOOZ  are  considered  daytime  data  and  are  plotted  in  red,  while 
data  collected  between  0100-1400Z  are  considered  niglittime  data  and  aie  plotted  in  blue. 
The  ah'  temperatme  collected  when  wind  speed  was  less  than  2  ms"^  are  also  ended  in 
gr  een.  It  is  clear  from  Figrne  33b  that  the  majority  of  the  static  tests  were  made  in  low 
wind  conditions.  Finthermore  Figrne  33b  shows  that  the  majority  of  the  overestimated 
temperatme  measinements  occiured  dining  low  wind  conditions. 

Static  T»st  Static  Test 


Static  Test 


•  Day 

•  Night 

®  WS<2ms‘’ 


Figrne  33.  Comparison  of  10  minute  averaged  static  Aumar  and  NFS  Met 
measinements  for  a)  pressine;  b)  air  temperatme;  c)  wind  speed. 
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Table  8  summarizes  the  statistical  comparison  between  the  Airmar  and  the  NPS 
Met  data.  Overall,  the  Airmar  follows  the  mean  trends  well,  but  may  deviate  considerably 
on  specific  instances.  The  Airmar  temperature  measurements,  in  particular,  appear  to 
suffer  from  solar  heating  of  the  sensor  housing  particularly  in  low  wind  conditions.  More 
testing  is  needed  to  better  understand  the  Airmar  performance  for  pressure  and  quantify 
the  temperature  error  associated  with  sensor  heating. 


Airmar  vs  NPS  Met  Static  Testing  (IHz  Data) 

Press  (mb) 

Temp  ( C ) 

Daytime  Temp  (  C  ) 

Nighttime  Temp  (  C  ) 

Wind  Sp  (ms"') 

Mean  Error 

2  94 

1  58 

2  13 

1  03 

0  17 

Error  Std  Dev 

0  58 

1  44 

1  78 

0  62 

0  39 

Abs  Error 

2  94 

1  76 

2  47 

1  04 

0  26 

Table  8.  Statistical  comparison  of  Airmar  and  NPS  Met  measurements 
using  10  minute  averages  of  the  1  Hz  raw  data  during  November 

2014  static  testing. 


C.  EVALUATION  OF  AIRMAR  MEASUREMENTS  COLLECTED  BY  LRI 

AND  NAVO  SHARCS 

The  Airmar  PB200,  as  the  default  METOC  sensor  on  the  LRI  METOC  Plus 
SHARCs,  has  been  used  in  multiple  at-sea  measurements  and  demonstrations.  Two  of 
these  data  collection  efforts  are  used  here  to  further  evaluate  the  Airmar  data  quality 
using  nearby  buoys.  The  LRI  PacX  SHARC  orbited  NDBC  buoy  46092  in  Monterey  Bay 
during  November  2011.  Two  NAVO  SHARCs  also  orbited  near  this  same  buoy  from 
September-December  2012.  The  proximity  of  SHARCs  to  the  NDBC  buoy  permits 
further  evaluation  of  the  SHARC  Airmar  PB200  at-sea  performance  using  the  hourly 
buoy  data.  This  data  was  interpolated  onto  the  time  of  the  SHARC  data  for  direct 
comparison.  The  results  are  provided  below. 

I.  PacX  SHARC  Data 

LRI  launched  four  SHARCs  from  San  Lrancisco  in  November  2011  to  begin  the 
PacX  event.  During  the  last  week  of  November  2011,  three  of  the  four  SHARCs 
remained  within  15  km  of  buoy  46092  for  approximately  four  days.  The  SHARCs  took 
turns  orbiting  the  buoy  and  then  loitering  nearby.  The  SHARC  sensors  should  be 

calibrated  at  its  best  at  the  time  of  this  data  collection  since  the  SHARCs  had  only 
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recently  departed  San  Francisco.  Figure  34  shows  the  LRI  PacX  SHARC  tracks  during 
this  timeframe.  The  inner  box  shows  the  tracks  that  are  within  15  km  of  buoy  46092. 


LRI  PacX  SHARCs  Nov  2011 
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Figure  34.  SHARCs  1,  2,  and  4  path  during  November,  2012  PacX.  The  inner 
box  shows  tracks  within  15  km  of  buoy  46092. 

Figure  35  is  the  time  series  plot  of  the  LRI  SHARC  10  minute  averaged  pressure 
and  the  NDBC  46092  hourly  averaged  pressure.  There  are  minor  gaps  in  both  the  NDBC 
buoy  data  and  the  SHARC  Airmar  data.  As  shown  in  the  NPS  Airmar  datasets,  the 
PB200  measurements  depict  the  pressure  trend  well  but  with  outliers  that  appear 
unreasonable.  Figure  36  shows  the  scatter  plot  of  the  SHARC  and  buoy  pressure  data. 
This  figure  also  reveals  that  the  two  measurements  are  generally  correlated  but  with  some 
“bad”  values  from  the  SHARC  Airmars. 
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LRI  PacX  /  Buoy  Pressure  Comparison 
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Figure  35.  Time  variation  of  the  PaeX  SHARC  and  buoy  46092  pressure 

measurements. 
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Figure  36. 


Buoy  46092  vs  LRI  PacX  SHARC  Measurements 


Comparison  of  the  LRI  SHARC  and  NDBC  pressure  measurements. 


Figure  37  shows  the  time  series  of  the  LRI  SHARC  and  NDBC  buoy  wind  speed 
and  wind  direetion  measurements.  The  buoy  measured  wind  speed  appears  to  undergo  a 
diurnal  cyele  that  the  Airmar  also  shows.  The  three  SHARCs  all  overestimated  wind 
speed  on  November  29,  2011.  The  reason  for  this  systematie  departure  is  not  known.  The 
wind  direetion  shifts  from  Southeast  to  North  during  this  PaeX  period  that  the  Airmar 
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largely  detects,  as  well.  Tlie  scatter  plots  (Figiue  38)  reveal  less  conelation  iii  these  data 
sets,  however,  than  in  the  NPS  at-sea  datasets  where  the  sensors  are  exactly  collocated. 


LRI  PacX  /  Buoy  Wind  Speod  Comparison 


Tims(UTC) 


Figure  37.  Same  as  Figiue  35,  except  for  a)  wind  speed;  b)  wind  dnection 

measiuements. 
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Buoy  46092  vs  LSI  PacX  SHARC  Measurements  guoy  46092  vs  LRI  PacX  SHARC  Measurements 


Figure  38.  Same  as  Figiue  36,  except  for  a)  wmd  speed  measiuemeuts;  b)  wind 

direction  measurements. 


The  buoy  ah  temperatiue  measiuemeuts  dming  this  time  period  vary  diiunally 
(Figiue  39).  The  SHARC  data  also  follow  this  pattern  but  appear  to  overestimate  the  peak 
temperatiues,  and  imderestimates  the  lows.  This  figiue  also  indicates  that  there  may  be 
multiple  factors  affecting  the  quality  of  the  Aumar  measiuement  resulting  in  the  sporadic 
consistency  shown  m  the  time  series  plot.  This  is  reflected  in  the  scatter  plot  as  well  with 
small  coiielations  of  the  two  measiuemeuts. 


LRt  PacX/  Buoy  Temperature  Comparison 
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Figure  39.  Same  as  Figiue  35,  except  for  au  temperatiue. 
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Buoy  46092  vs  LRI  PacX  SHARC  Measurements 


Figure  40.  Same  as  Figure  36,  except  for  air  temperature  measurements. 

To  further  analyze  the  SHARC  Airmar  temperature  measurements  and  examine 
the  effect  of  solar  heating  on  the  resultant  data,  I  separated  the  data  into  daytime  and 
nighttime  observations  following  the  same  method  used  on  the  NFS  static  data.  Again, 
the  air  temperature  that  coincided  with  wind  speed  less  than  or  equal  to  2  ms'*  was  circled 
in  green.  As  seen  in  Figure  41a,  the  vast  majority  of  the  overestimated  Airmar 
temperatures  occur  during  the  day.  Figure  41b  reveals  a  direct  association  between 
overestimated  air  temperature  data,  daytime  heating,  and  low  wind  speed.  This  plot 
indicates  that  the  solar  heating  effect  on  the  Airmar  can  be  significant. 
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Buoy  46092  i/s  LRI  PacX  SHARC  Measurffnents 


Buoy  46092  vs  LRt  PacX  SHARC  Measurements 


Figure  4 1 .  Comparison  of  the  LRI  PacX  SHARC  and  NDBC  buoy  46092 

measmed  an  temperatme  with  a)  daytime  data  in  red  and  nighttime  data 
in  blue;  b)  when  wind  speed  was  less  than  2  ms'*  (circled  in  green)  to 
examine  the  effects  of  solar  heating  on  the  Airmar. 


2.  NAVO  SHARC  Airmar  Data  Analysis 

From  September-December  2012,  NAVO  steered  two  of  thens  SHARCs  to  orbit 
in  Monterey  Bay.  During  this  tune,  the  NAVO  SHARCs  were  within  15  km  of  the 
NDBC  biroy  46092.  The  NAVO  SHARC  tracks  are  showrr  in  Figine  42.  NAVO  powered 
down  the  SHARC  AitTuar  weather  stations  for  long  periods  of  time  to  conserve  power  (S. 
Crossland  2014,  NAVO,  personal  commimication)  birt  did  collect  approximately  12  days 
of  near  sirrface  measinements.  Figine  43  shows  the  wind  direction  from  the  entire  period 
to  illustrate  the  presence  of  the  Airmar'  weather  data  gaps.  Of  note,  the  buoy  wind 
direction  resolution  is  10  degrees  while  the  SHARC  Airmar  resolution  is  1  degree.  The 
scatter  plot,  Figine  44,  reveals  large  number  of  correlated  measmements  from  the  two 
systems,  although  we  can  see  clearly  a  large  number  of  imconelated  results  as  well. 
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Figure  42.  NAVO  SHARC  tracks  within  15  km  of  buoy  46092  from  September 
-  December  2012  in  Monterey  Bay.  Location  of  Buoy  46092  is  shown 
in  green  on  the  horizontal  axis. 


NAVO  SHARC  /  Buoy  Wind  Dir.  Comparison 


Figure  43.  Time  variation  plot  of  the  NAVO  SHARC  Airmar  and  NDBC  buoy 
46092  wind  directions  from  the  entire  intercomparison  mission  between 
September  and  December  2012. 
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Buoy  46092  vs  NAVO  SHARC  Measurements 


Figure  44.  Comparison  of  the  NAVO  SHARC  and  NDBC  buoy  46092  wind 

direction  measurements. 

For  the  next  set  of  time  series  plots,  the  time  periods  in  October  and  November 
2012  when  both  the  NAVO  SHARC  and  NDBC  buoy  were  collecting  data  are  isolated  to 
show  more  details  (Figure  45).  Figures  45a  and  45b  contain  the  wind  speed  data  collected 
by  the  NAVO  SHARCs  and  NDBC  buoy  46092  in  October  and  November,  respectively. 
Figures  45c  and  45d  show  the  pressure  measurements  for  these  platforms  in  October  and 
November,  respectively.  Both  the  sets  of  time  series  plots  and  the  scatter  plots  shown  in 
Figures  46a  and  46b  indicate  that  the  SHARC  Airmar  data  have  generally  good 
correlation  with  the  NDBC  buoy  data.  There  are,  however,  outliers  and  the  SHARC 
pressure  measurements,  in  particular,  are  prone  to  have  outliers.  Again,  these  plots 
indicate  that  the  Airmar  provides  good  trend  information  but  not  the  details,  and  require 
data  quality  control. 
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NAVO  SHARC I  Buoy  Wind  Sp«»d  Comparison  NAVO  SHARC  t  Buoy  Pressure  Comparison 


Tim*  (Day.  UTCi  Time  (Day.  UTC) 


Figure  45.  Time  variation  of  the  NAVO  SHARC  and  NDBC  buoy  46092 

measurements  with  a)  wind  speed  centered  aroimd  October  1,  2012;  b) 
wind  speed  centered  aroimd  November  1,  2012;  c)  pressiue  centered 
around  October  1,  2012;  d)  pressme  centered  around  November  1, 

2012. 


Buoy  4«092  vs  NAVO  SHARC  Messurcmsnts  Buoy  46092  vt  NAVO  SHARC  M«aiur«mtn(s 


Figure  46.  Compar  ison  of  the  NAVO  SHARC  and  NDBC  buoy  46092  a)  wind 
speed  measiuements;  b)  pressme  measurements. 
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Figiue  47  shows  the  time  series  plots  of  the  NAVO  SHARC  and  NDBC  buoy 
46092  ah  temperatiue  measiuemeuts  centered  about  October  1,  2012  and  November  1, 
2012,  respectively.  The  NAVO  SELA.RC  data  follows  the  general  dimiial  cycle  but, 
overall,  the  SHARC  overestimates  the  ah  temperatiue.  Following  the  same  process  used 
with  the  LRI  PacX  data  to  examine  the  effect  of  solar  heathig  on  the  Airmar,  I  separated 
the  day  and  night  observations  (Figure  48a)  and  then  chcled  data  that  corresponded  to 
periods  of  low  wind  speed  (Figiue  48b).  While  the  overly  high  temperatiues  are  not  as 
well  correlated  with  daytime  heathig  errors,  they  are  strongly  associated  with  low  wind 
speed. 


Figirre  47.  Same  as  Figiue  45,  except  for  ah  temperatme. 


Figiue  48.  Same  as  Figiue  41,  except  for  the  NAVO  SHARC. 
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3. 


Airmar  Data  Set  Conclusions  and  Statistics 


Overall,  the  LRI  PaeX  and  NAVO  Airmar  datasets  aeeurately  depiet  longer  seale 
trends.  The  data  requires  quality  eontrol,  however.  The  Airmar  temperature  sensors  seem 
to  have  different  bias  day  and  night,  suggesting  the  role  of  solar  heating  resulting  in  the 
additional  errors.  Table  9  summarizes  the  statistieal  evaluation  of  Airmar  measurements 
against  other  known  sensor  paekages.  Of  note,  the  LRI  and  NAVO  Airmar  have 
temperature  mean  errors  of  2.02  and  1.00  °C  during  eonditions  of  daytime  heating  and 
low  winds.  To  eombat  these  issues,  the  Airmar  sensor  should  be  ealibrated  separately  for 
daytime  and  nighttime.  Speeifieally  during  the  daytime,  the  Airmar  should  be  ealibrated 
against  a  high  quality,  temperature  sensor  in  aspirated  radiation  shield.  Sinee  the  Airmar 
does  not  have  a  well-ventilated  radiation  shield,  the  magnitude  of  solar  heating  effeet  is 
greater  than  that  experieneed  by  the  NPS  Met  and  NDBC  buoy  sensors. 
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Airtnar  vs  NFS  Met  At  Sea  (10  Minute  Averages) 

Press  (nib) 

Temp  (  C  ) 

Wind  Sp  (ms’’) 

Mean  Error 

-0.39 

3.28 

0.09 

Error  Std  Dev 

0.21 

1.00 

0.19 

Abs  Error 

0.40 

3.28 

0.18 

Corr  Coef 

0.81 

0.64 

1.00 

NAVO  SHARCs  vs  NDBC  Buoy  46092 

Press  (n±)) 

Wind  Sp  (ms  *) 

Day  Temp  witliWmds  <2  (ins  ’) 

Mean  Error 

1.08 

0.36 

2.02 

Error  Std  Dev 

1.47 

1.41 

0.88 

Abs  Error 

1.50 

1.13 

2.03 

Corr  Coef 

0.89 

0.82 

0.65 

Temp  (  C ) 

Daytime  Tenp  (  C  ) 

Ni.dittime  Temp  (  C  ) 

Mean  Error 

0.75 

0.84 

0.68 

Error  Std  Dev 

1.06 

1.20 

0.93 

Abs  Error 

1.05 

1.16 

0.96 

Corr  Coef 

0.75 

0.70 

0.80 

LRI  PacX  SHARCs  vs  NDBC  Buoy  46092 

Press  (n±)) 

Wind  Sp  (ins  ’) 

Day  Temp  witliWmds  <2  (ms  ’) 

Mean  Error 

0.34 

1.07 

1.00 

Error  Std  Dev 

1.43 

2.00 

1.31 

Abs  Error 

0.73 

1.61 

1.31 

Corr  Coef 

0.96 

0.68 

0.80 

Temp  (  C ) 

Daytime  Tenp  (  C  ) 

Ni.dittime  Temp  (  C  ) 

Mean  Error 

-0.57 

-0.32 

-0.97 

Error  Std  Dev 

1.33 

1.51 

0.83 

Abs  Error 

1.18 

1.27 

1.04 

Corr  Coef 

0.73 

0.72 

0.81 

Table  9.  Statistical  comparison  between  SHARC/Airmar  at-sea 
environmental  measurements  and  measurements  from  known,  high 

quality  instruments. 
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V.  EXAMPLE  OF  SHARC-BASED  NEAR  SURFACE 
MEASUREMENT  FOR  NAVAL  APPLICATIONS 


In  addition  to  validating  the  quality  of  the  NPS  Met  sensors,  the  suitability  of 
using  the  NPS  Met  near  surface  mean  observations  for  Naval  applications  is  examined 
here.  The  NPS  Met  at-sea  measurements  were  used  as  inputs  to  the  COARE  algorithm  to 
calculate  surface  fluxes  of  momentum,  sensible,  and  latent  heat,  the  near  surface  mean 
thermodynamic  profiles,  and  the  EDH.  This  step  also  further  evaluates  the  suitability  of 
the  SHARC  as  a  platform  for  obtaining  high  quality  measurements  for  input  into  coupled 
forecast  models  and  EM  propagation  prediction  models. 

A.  SURFACE  FLUX 

THE  NPS  Met  mean  air  temperature,  SST,  wind  speed,  and  RH  were  input  into 
the  COARE  algorithm  to  derive  parameterized  momentum,  sensible  heat,  and  latent  heat 
fluxes.  These  derived  fluxes  were  then  compared  to  the  turbulent  fluxes  measured  by  the 
MASFlux  buoy.  The  MASFlux  buoy  includes  a  sonic  anemometer,  mounted  at  3.48  m, 
which  directly  measures  virtual  temperature  and  3 -dimensional  wind  components  at  a 
sampling  rate  of  20  Hz,  which  yield  virtual  sensible  heat  and  momentum  fluxes  using  the 
eddy  correlation  (EC)  method. 

Figure  49  shows  the  time  series  plots  of  NPS  Met  and  MASFlux  surface  fluxes 
derived  using  mean  measurements  and  COARE,  and  the  fluxes  from  the  EC  method,  for 
each  day  of  at-sea  testing.  The  derived  fluxes  from  the  NPS  Met  and  MASFlux  compare 
well  to  each  other  and  any  differences  in  flux  are  due  solely  to  differences  in  the  mean 
quantities.  The  derived  fluxes  show  the  same  generally  increasing  trend  as  the  EC  fluxes 
on  July  29  and  July  30,  2014,  but  do  not  capture  as  much  variability  as  the  directly 
measured  fluxes  (Figures  49a  and  49b).  On  August  8,  2014,  the  MASFlux  3-D  sonic 
measured  momentum  fluxes  are  almost  50  percent  higher  than  the  fluxes  derived  from 
the  NPS  Met  and  MASFlux  mean  measurements  (Figure  49c).  This  difference  is  likely 
due  to  the  high  swell  conditions  encountered  that  day.  The  buoy,  the  SHARC,  or  both 
platforms  may  have  been  in  the  wave  boundary  layers  where  basic  assumptions  of  MOST 
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and  COARE  aie  violated.  If  the  platfomis  were  within  the  wave  boundary  layers,  this 
factor  would  also  negate  the  assumption  of  constant  fluxes  within  the  siuface  layer. 


20140729  SHARC  vs  MASFlux 
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Figure  49.  Time  variation  of  NFS  Met  and  MASFlux  siu  face  fluxes  derived 
using  mean  measiuements  and  COARE,  and  the  fluxes  from  the  EC 
method  on  a)  July  29,  2014;  b)  July  30,  2014;  c)  August  8,  2014. 


Figiues  50  is  the  scatter  plot  of  NFS  Met  parameterized  momentum  flux  versus 
that  from  the  MASFlux  using  the  EC  method.  Similarly,  Figiues  5 1  is  the  scatter  plot  of 
NFS  Met  parameterized  vutual  sensible  heat  flux  versus  that  from  the  MASFlux  using 
the  EC  method.  Although  the  mean  measiuements  and  COARE  tended  to  underestimate 
the  fluxes,  both  of  these  plots  show  generally  good  conelation  between  the  derived  and 
measiued  fluxes  on  July  29,  2014  and  July  30,  2014.  The  conelation  is  poor,  however,  on 
August  8,  2014.  Again,  if  the  SHARC  platform  was  in  the  wave  boundary  layers,  then 
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COARE  would  be  insufficient  for  deriving  flux  values  from  mean  measurements  and 
hence  the  poor  correlation  with  EC  fluxes  shown  in  Figures  50  and  51. 


MASFlUX  EC  VS  SHARC  COARE 


Figure  50.  Comparison  of  NFS  Met  momentum  flux  derived  using  mean 
measurements  and  COARE,  and  the  momentum  flux  from  the  EC 

method. 


MASFlUX  EC  vs  SHARC  COARE 


Figure  5 1 .  Same  as  Figure  50,  except  for  virtual  sensible  heat  flux. 
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B,  NEAR  SURFACE  PROFILES  AND  EDH  CALCULATIONS 

In  addition  to  surface  fluxes,  the  SHARC  mean  measurements  and  the  output  of 
the  COARE  algorithm  ean  be  used  to  derive  the  surfaee  layer  vertieal  profiles  of  mean 
wind,  potential  temperature  and  speeifie  humidity.  The  thermodynamie  profiles  ean  then 
be  used  to  determine  the  vertieal  profile  of  the  refraetivity  index  or  the  modified 
refraetive  index,  M  (Eq.  9),  which  determines  EM  ducting  and  propagation  paths  in  the 
surfaee  layer.  Figure  52  shows  an  example  of  two  profiles  caleulated  using  the  mean 
measurements  eollected  by  the  SHARC  and  NPS  Met  sensors  on  August  8,  2014.  As 
deseribed  in  Chapter  II,  the  EDH  is  one  measure  of  the  property  of  the  evaporative  duct 
and  its  ability  to  trap  EM  waves  and  is  defined  by  the  height  at  whieh  dM/dz  is  equal  to 
zero.  The  evaporation  duet  strength  (EDS)  is  defined  as  the  M-  differenee  between  the 
EDH  level  and  the  surfaee.  Figure  53  shows  a  plot  of  the  EDH  and  EDS  at  10  minute 
intervals  on  August  8,  2014  as  determined  using  SHARC  mean  measurements  and  the 
derived  M  vertical  profile.  The  red  and  blue  lines  indieate  the  time  when  the  example 
profiles  in  Figure  52  are  taken  and  the  assoeiated  EDH  and  EDS  at  these  time  instanees. 
The  SHARC ’s  ability  to  colleet  quality  near  surfaee  measurements  allows  users  to  obtain 
near  real  time  vertieal  profiles  and  the  assoeiated  EDH  and  EDS  quantities  which  can  be 
input  to  the  atmospherie  propagation  models,  sueh  as  AREPS,  for  predietion  of  EM 
propagation  as  a  taetieal  deeision  aid. 
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0  q  M 

Figure  52.  Example  vertical  profiles  derived  using  SHARC  measurements  from 
August  8,  2014  at-sea  testing  and  the  output  of  COARE  algorithm. 


77 


Time  (Hour  UTC) 


Figure  53.  EDH  and  EDS  calculated  from  the  SHARC  derived  near  surface 
vertical  profdes  during  August  8,  2014  at-sea  testing. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SHARC,  NPS  MET,  AND  AIRMAR  CONCLUSIONS 

The  SHARC  is  a  capable,  suitable,  and  advantageous  moving  platform  for 
collecting  near  surface  measurements.  The  SHARC  platform  offers  numerous  new  and 
exciting  capabilities,  including  its  high  endurance  and  ability  to  withstand  harsh 
conditions.  Additionally,  operators  can  steer  the  SHARC  to  areas  of  interest  and  have  it 
remain  on  station.  Critically,  this  thesis  also  shows  that  the  SHARC  platform  does  not 
distort  the  mean  measurements  collected  from  it. 

The  NPS  Met  is  proven  to  be  better  suited  for  obtaining  high  quality  mean 
measurements  for  air-sea  interaction  studies  and  Naval  applications  than  the  default 
Airmar  PB200  sensor  package.  The  NPS  Met  and  the  known,  high  quality  MASFlux 
mean  measurements  are  highly  correlated  and  display  only  small  biases  due  to  spatial  and 
temporal  variability.  The  NPS  Met  package  also  measures  RH  and  SST  which  are  not 
available  via  the  Airmar,  and  are  critical  to  determining  surface  fluxes,  near  surface 
profiles,  and  properties  of  the  EDH.  The  low  power  consumption  of  the  NPS  Met  also 
lends  itself  well  to  long  endurance  missions  on  the  SHARC.  Although  solar  radiation 
heating  of  the  air  temperature  sensor  does  affect  the  NPS  Met  and  all  naturally-ventilated 
systems,  the  radiation  shield  minimizes  this  effect. 

The  Airmar  PB200  captures  mean  trends  adequately  but  requires  individual 
instrument  calibration  and  data  quality  control.  The  Airmar  temperature  sensor  requires 
careful  calibration,  both  to  eliminate  biases  and  to  account  for  solar  heating  errors.  As  the 
Airmar  does  not  have  an  effective  radiation  shield,  the  magnitude  of  solar  heating  is 
greater  than  that  experienced  by  the  NPS  Met  and  NDBC  buoy  sensors,  and  the  effect  is 
especially  evident  in  low  wind  conditions.  The  Airmar  pressure  measurements  appear  to 
be  inconsistent  at  times  and  also  require  data  quality  control.  Finally,  the  Airmar  does  not 
sample  RH  or  SST  which  are  required  for  many  Naval  applications. 
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B, 


RECOMMENDATIONS  FOR  FUTURE  WORK 


In  order  to  maximize  SHARC  and  NPS  Met  usages  and  functioning,  several  areas 
of  research  require  more  work.  For  groups  continuing  to  use  the  default  Airmar  weather 
instrument,  there  needs  to  be  a  more  quantitative  process  or  algorithm  to  provide  pressure 
and  temperature  data  quality  control.  Additionally,  the  suitability  of  the  SHARC  for 
measuring  waves  and  sea  state  should  be  evaluated.  As  part  of  this  process,  the  wave 
sensor  recommended  by  LRI  SHARC ’s  METOC  Plus  model  should  be  further  evaluated. 
The  NPS  Met  and  MASFlux  measurements  and  data  sets  available  from  NAVO  provide 
an  excellent  starting  point  for  this  task.  For  the  NPS  Met  package  and  SHARC  platform, 
future  research  should  focus  on  the  effect  of  the  wave  boundary  layer  on  mean  and  flux 
measurements,  especially  during  high  swell  or  wave  conditions.  Finally,  the  NPS  team  is 
in  the  process  of  developing  a  direct  flux  measuring  payload  for  the  SHARC.  This  new 
payload  will  require  testing  and  evaluation,  as  well. 
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